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Abstract: 

Q^ | In honor of Onsager's ninetieth birthday, we like to review some exact 

results obtained so far in the chiral Potts models and to translate these results 
into language more transparent to physicists, so that experts in Monte Carlo 
calculations, high and low temperature expansions, and various other 
methods, can use them. 

We shall pay special attention to the interfacial tension e r between the k 

Cn ' state and the k — r state. By examining the ground states, it is seen that the 

integrable line ends at a superwetting point, on which the relation e r = re\ 
is satisfied, so that it is energetically neutral to have one interface or more. 

\^0 ■ We present also some partial results on the meaning of the integrable line for 

low temperatures where it lives in the non-wet regime. We make Baxter's 
exact results more explicit for the symmetric case. By performing a Bethc 
Ansatz calculation with open boundary conditions we confirm a dilogarithm 
identity for the low-temperature expansion which may be new. 

q«^ | We propose a new model for numerical studies. This model has only two 

variables and exhibits commensurate and incommensurate phase transitions 
and wetting transitions near zero temperature. It appears to be not 
integrable, except at one point, and at each temperature there is a point, 
where it is almost identical with the integrable chiral Potts model. 

a: 

KEY WORDS: Chiral Potts model; chiral clock model; star-triangle 
equations; Yang-Baxter equations; interfacial tension; wetting; superwetting; 
scaling; corrections to scaling; low-temperature expansions; dilogarithms; 
Bethe Ansatz. 
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1. INTRODUCTION 

When Onsager published his solution of the two-dimensional Ising model 
in 1944,® this was almost instantly recognized as a milestone in the 
development of statistical mechanics. Many new developments were inspired 
by his results. On the other hand, Onsager's techniques were far ahead 
of his time and, when he announced his incredibly simple result for the 
spontaneous magnetization as a comment to a conference talkpcP his paper 
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attained a mystical status for many years after. That his student Kaufman 
simplified the solution through Clifford algebras (i.e. fermions)®tP did little 
to change this. 

Now, nearly fifty years later, we can appreciate Onsager's methods 
much better. He was the first to introduce the star-triangle equation® to 
statistical mechanics^® even though now the name Yang-Baxter equation is 
commonly used.CrE® He also was the first to introduce loop algebras O as a 
solvability principle. Onsager and Kaufman have clear priority over Wick for 
the Wick theorem® and, together with Bethe,^ thev scouted a new area of 
mathematics which is now called quantum groups. (£L34LZl> 

Onsager's (only partly published) work on two-point functions in the Ising 
model got extended in 1966,t2r El but it was only in 1973, when Wu, McCoy, 
Tracy, and Barouch^ErEif announced the Painleve equation for the scaled 
two-point correlation, that the theory of the two-dimensional Ising model 
went beyond Onsager's level. Also, the first two-dimensional models solved 
that were more complicated than the Ising model were Lieb's ice model^EAE^ 
of 1967 and Baxter's eight-vertex model® of 1973. 

Onsager's loop-algebra solution method was generalized only in 1985 
when Von Gehlen and Rittenberg^P solved the Dolan- Grady® criterium 
within a one-dimensional generalization of the quantum Pottsmodel. The 
connection with Onsager's 1944_paper was noticed by Perk/Bat^P showing 
that the chiral Potts modelfflEl'El) is the first genuine generalization of the 
Ising model, with its "superintegrable" case® solvable for two reasons: star- 
triangle integrabilitySEllilf and loop-algebra integrability.Sc^P The chiral 
Potts model upgrades the fermions of Kaufman® to parafermionsll It also 
provides an infinite hierarchy of quantum groups at roots of unity, with Ising 
as its first entry®© Thus the works of Onsager are still at the center of 
attention. 

This paper is organized as follows: In section |2| we define the chiral Potts 
model and review some of the exact results obtained so far.^EZI - ^ We pay 
special attention to the interfacial tensions e r , giving several results that 
have not been presented before as such.a We show that the solvable chiral 

3 We note that the parafermions here are of "cyclic root-of-unity" type, 1 © generalizing 
the Weyl algebra, not of the more usual "highest- weight" type introduced by Green. v3 
Within the original Ising model and its fermion approaches these two types are isomorphic, 
however. 

4 We gratefully acknowledge several private communications with Dr. Baxter. 



Potts models "superwet," that is e r = re±, at T = and T = T c . However, 
BaxterHH has shown that, for < T < T c , the interfacial tensions satisfy 
the inequalities e r < e r _j + €j, so that the integrable line is in the not-wetted 
region. We use low-temperature expansions and the exact results of Baxter 
to further analyze the effects of various boundary conditions. We also discuss 
what the critical exponents obtained exactly mean for the scaling function. 
In sections [| to |5], we discuss the relation between the integrable subcases, 
whose Boltzmann weights can be represented by a product-form, and the 
generalized clock-model representations of the pair-interaction energies. 
Special attention is given to the symmetric case (with equal horizontal and 
vertical interactions) and several numerical details and graphs are given. In 
section || we outline a Bethe Ansatz calculation for the leading corrections 
to the zero-temperature diagonal interfacial tensions. Finally, in section [7|, 
we introduce a new model with only few parameters, which is very close to 
the integrable model and may deserve detailed further study by numerical 
means. 

2. CHIRAL POTTS MODEL 

2.1. Integrable Chiral Potts Model 

In our original paper,® new exact-solution manifolds were discovered within 
the chiral Potts model. To be more specific, the Potts model,^2rE§ which is 
itself a generalization of the two-dimensional Ising model solved by Onsager w 
in 1944, and whose interaction energy for the two spins on an edge is given 

by 

£(a,a')=E5 a , a ,, (2.1) 

was generalized to the chiral Potts model (or Z g -modelteP) with interaction 
energy 

N-l 

£{cr, a') = £(n - n') = £ % uj j{n - n '\ (2.2) 

3=1 

where 

, , 2m/N „ , ,n / , ,n' /n q\ 

oj = e ' , a = uj , a — u) . {^■^) 



Clearly the interaction energy defined by ( |2.2| ) satisfies the relation S(a, a') 
£{n - n') = 8{n - n' + N). For Ej = E, (1 < j < N - 1), using 



N-l 



J- u,i(— ') = NS ff y - 1, (2.4) 



=i 



we find that the interaction energy is identical to that of the Potts model 
except for an overall constant. When Ej = E^-j, it includes the integrable 
Fateev-Zamolodchikov^) self-dual Z^r model as a special case. 
The Boltzmann weight for an edge is 

W(n - ri) = e~ £{n ~ n ' )/k * T (2.5) 

with &b the Boltzmann constant. Since duality transform is equivalent 
to Fourier transform^E^ of these weights, the weights are said to be 
"self-dual," if they are equal for Droportional) to their Fourier transforms. 
With a great deal of effort Hc§EP and some luck we were able to find 
self-dual solutions to the star-triangle equations and found that they can 
be written in product forms. 

In Australia, Baxter and the two of usW - - mainly through guessing 
and especially guided by Onsager's work^ -- found more general solutions 
of the star-triangle or "checkerboard Yang-Baxter" equation, 

N 

Y,W qr (b-d)W pr (a-d)W pg (d-c)=R pqr W pq (a-b)W pr (b-c)W qr (a-c). 

d=l 

(2.6) 
They are also given in product form. To be more specific, we found that 
in the chiral Potts model, the weights W pq or W pq depend on two line (or 
rapidity) variables denoted by p = (x p ,y p ,fi p ) and q = (x q ,y q , fi q ), shown in 
Fig. 1; the weights are given as 

W pq {n) _ ffi p \ (y q -x p u n \ W pq {n) _ ,f ujx p 



W pq (n - 1) V/V \y P - x q uj n )' W pq (n - 1) ' "' ' V y q - y p uj n 

(2.7) 
Here, the parameters p and q are restricted by the two periodicity 
requirements W pq (N + n) = W pq {n) and W pq (N + n) = W pq (n), yielding 
the conditions 

/,, \N y N - X N V N -V N 




q 




Fig. 1. Boltzmann weights W pq (a — b) and W pq (a — b) for two types of edges connecting 
spins a and b. As W(a — b) ^ W(b — a), and similarly for W, we need to put arrows on 
the edges to distinguish the two different choices. 



These imply the existence of numbers k and k' related by A; + k' 
that the equations 



1 such 



/ J 



N 



k'/(l - kx N p ) = (1 - ky?)/k' 



*Z+v2 



Ml + «) (2-9) 



hold for each p. (This may require using the ambiguity in defining the 



x 



v- 



x q , y p , and y q in 



and rescaling them with a common factor.) The 
equations (|2.9|) describe a complex curve, and the genus of this curve is 
g = N 2 (N -2) + l. 

From Fig. 2, one can see that the star-triangle equations allow one to 
move the rapidity line p through the vertex (the intersection of the other 
two rapidity lines). Because of this, one can permute these rapidity lines 
without changing the partition function, except possibly some constant 
factors. Baxter called such lattice models Z-invariant.© This also means 
that transfer matrices associated with different rapidity variables commute. 
We can see from (|2.7p to (|2.9|) that for given k, there is only one free 



variable x p , associated with each rapidity variable p, and that y p and fx p 



can be determined from ( |2.9| ). For the rectangular lattice with just two 
rapidity variables p and q, and therefore two kinds of weights W pn and W r . 
there are three free variables. 
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Fig. 2. The star-triangle relations, allowing one to move the rapidity line p through a 
vertex, which is the intersection of the two other rapidity lines q and r. 

By comparing with the Ising model (JV = 2),® where k and k! are the 
elliptic modulus and its complementary modulus, we conclude that k and k' 
describe how far the system is from its critical point. 

2.2. Gauge Transformations 

Moreover, let 



Vp — V[ x p> "v> 



W pq (n) 



W pq {n), 



W vq {n) 



(VpVqT W pq (n) 



(2.10) 



for any arbitrary function r\. We can then replace the W and W in (|2.6|) by 
the W and W' inJ^Tq). We find that, if W pq (n) and W m {n) satisfy Q, 
then W pq {n)' and W pq {n)' also satisfy the star-triangle equations ( |2.6| ). 

Furthermore, the transformation ( p.lOj ) leaves the partition function for a 
system with periodic boundary conditions invariant. This can be seen easily 
by examining what happens at a particular site e under such a transformation 
in a checkerboard lattice with p, p', q, and q' as the rapidity variables, as 
shown in Fig. 3. One finds that the additional factor if in W pq (e — d)' cancels 
out the factor r]~ e in W p i q [e — c)', and the rf in W pq (e — d)' cancels out the 
r]~ e in W pq '(a — e)', etc, leaving the net contribution at each site unchanged. 
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Fig. 3. Gauge transformation at a particular site e adds gauge factors to each of the four 
Boltzmann weights associated with the four bonds meeting in e. The total contributions 
add up to zero. 



then W pq (n) 



and W pq (n)' are no 



Particularly, if we choose r\ p = fi~ 
longer periodic, and they differ from (|2.7|) by dropping the /i p and \x q factors, 
making the weights and transfer matrices depend rationally on x p , x q , y p , and 
y q only, and therefore more manageable. Indeed, as we are left only with the 
last equation in ( |2.9| ) the genus of the curve is then reduced to g = (N — l) 2 . 

On the other hand, if 7] (x, k) in ( p.lOj ) is a constant function, then only 
W picks up a factor. When we let r\ p = r\ q = u p , then W remains periodic, 
whenever p is an integer. 



2.3. Integrable Model and Zero-Temperature Limit 

We now compare the weights of the integrable model with the weights given 
by fl21]) and Q, which has 2(N- 1) variables Ej and Ej. Hence for (^) to 
be integrable, there must be 2N — 5 equations between these 2N — 2 variables. 
We rewrite the 2(N — 1) variables as 



E^ 
k B T 



Kjtu^, 



E 



N- 



&=*"*' 



k B T 

En-j 
~kB~T 



Kjuj- a >, 



Kjuj-^, 



(2.11) 



for 1 < j < [N/2], with [x] denoting the integral part of x. Therefore ( |2.2| ) 
becomes 



' |(JV-1) 



g(») 



/2-7T \ 

£ 2i^cos^— {jn + Aj)), Nodd, 



N-l 



£ 2Kj cosQijn + A,) ) + K hN {-iy 



N even, 



(2.12) 

with similar equations for 8. For real Kj and Aj, the Boltzmann weights 
are real and positive. When N is odd, we can think of the interactions as 
composed of ~(iV — 1) chiral clock model terms; particularly, for N = 3, 
it is the three-state chiral clock model. For even iV there is an additional 
Ising-like term; for example, for N — 4, it is composed of a four-state clock 
model with an Ising term. 

The weights of the integrable models given in ( |2.7| ) through ( |2.9| ) can be 
rewritten in the form 



W(n) ((!,«))„,„ W(n) ((l,a)) , 



W(0) ((l,/3)) ,n' 

where we have used the definitions 



W(0) ((l,/?)) ,n 



(2.13) 



1, (y.))ran 



A(a) n - m 



A(a) 



'l-a 



N\VN 



l,a) m ,„ = J| (l-uj^a), 

j=m+l 

(2.14) 



for m < n. It is easy to verify that 



1, QiJJ^f! 



'l,a) 



m,,fe' 



'l,a) 



fc,n- 



Since ((1, a)) m , m = 1, we can extend the definition to m > n by 

((1, a)) m ,„ = 1/((1, a))«, m . 
Moreover, because of the normalization factor A(a), we have 

((1, Oi)) mm+N = 1. 



(2.15) 



(2.16) 



(2.17) 



Hence, the weights in (|2.13|) are always periodic with period N. 



In order that the weights ( |2.13| ) satisfy the star-triangle equations 
they have to satisfy only one necessary condition on the four constants a, (3, 
a, and /3, namely 

~ R = ^- 2.18 

p to a 

This is consistent with ( |2.7|) ; it can also be derived directly. Substituting 
W{n), W(n) (for W pq and W pq ), W'(n), and W\n) (for W^. and W pr ), as 
given by ( |2.13| ), into ( |2.6| ), we can solve W"(n), and W"(n) (for W qr and 



Wqr), provided 



a 

7 







uj a 



a' 

J' 






A(a)A(a) A(aOA(aO 



A(/3)A(/3) A(J3')A(J3') 



(2.19) 



From this, we precisely reproduce the Z-invariant periodic solutions (|2.7| ) to 
( |2.9|) , up to possible gauge transformations as mentioned in (|2.10| ). Other 
than that, the only other allowed variation on the weights is that u may 
be chosen to be any root of unity tu N = 1 or w = e 2m j/N f QT an y j n ^ e g er 
1 < j < A^ — 1. Here j and A^ must be relative prime, otherwise we would 
need to redefine A(a) in order to retain periodicity. Letting 



a 



,.2'ie 



P 



2i(j> 



a: 



„2i8 



8 



D 2i0 



(2.20) 



then the weights can be rewritten as 



W{n) 
W(0) 

W(n) 

W(o) 



n/N 



sin(A^0) 
sm(N6) 



sin(iV0)^ /A 



sin(6» + nj/N) 



sm(N6)_ 




(2.21) 



which are real as long as 9, 



sin(0 + irj/N) 
9, and 4> are real. Now ( [2.1 8|) becomes 
9 — -(f) + 9. (2.22) 



7T 

N 



Therefore, the weights are functions of three independent variables. We can 
relate the "elliptic modulus" k with these three variables as most of the exact 
results are given in terms of this k. Comparing (|2.7|) with fl2.13|) , we find 



a 



„2i0 



'' i. 



Vi 



P 



D 2J0 



,L r , 



y P 



a = e 2ie 



J ( . 



UJXr 







2i<f> 



Vp 



(2.23) 



Therefore 




e 2i(6-4>) _ x pVp e 2i{4> 

•EqVq 


We can then use (|2.9|) to obtain 


k 2 = sin" 2 (N(6 - 0)) 


cos 2 



^ x p 2i{6+8) _ X q 



V P uy q 

(n(9 + 9)) + cos 2 (N(<f> - 9)) 
-2 cos (N(9 + 9)) cos (N(6 - 6)) cos (N(9 - 0))] . (2.25) 



This is also easily verified substituting (|2.24j ) in the right-hand side of ( [2.25| ) 
and eliminating y^ and y^ using the last equality in 



In section |3|, we express 9 and 6 i n ( P-20|) and ( |2.21| ) in terms of the N — 1 
variables Kj and Aj of ( |2 . 1 2| ) . Clearly, these variables must satisfy N — 3 
consistency relations. The four variables 9, 6 and 9, 6 can be easily rewritten 
in terms of Kj, Aj, Kj, Aj. The integrability condition (|2.22|) gives another 
condition relating them. 

For the square lattice with W = W, (which is called the symmetric case 
in the following), it then follows from ( [2.2 1| ) that we must have 9 = 9 and 
= 0. Consequently, the integrability condition becomes 

6-9=— or 6 = 9 + —. (2.26) 

v IN v 2N K J 

Now it is very easy to express Kj = Kj and Aj = Aj in terms of the single 
variable 9, and plot graphs for different N. The details are included in section 
|^ and here we outline a few of the conclusions. We find, with the energy unit 



convention k^TKi = 1 of subsection [572], that the integrable curve ends at 
zero temperature at 



AT-l 



sin(7r/iV) f2imj A j x 

M^jN) cos {— ±( ] 

sm(n/N) JnT(2l _ N)/4 



~ s{n) = g m^JJn) cos VN~ ± (2 " W 

N-l 

= y — — —u 

fr! sm(nl/N) 



i[ 



N±2n-l + 2N[ T §}) sin ^. (2.27) 

The last form is linear in n, periodically extended, with [x] again denoting 
the integral part of x. The first form in ( |2.27| ) hides the linearity with n, but 

10 



will suggest an interesting generalization. As the temperature T increases, 
Ai decreases, and on the integrable line, the ratios Kj/ K\ and Aj/Ai remain 
almost constant for 1 < j < [h(N — 1)], namely 

Kj _ sin^/AQ A,- _ N - 2j 

K[ ~ sin^/iV) + K » AT" A^2" + ^' (2 ' 28) 

with Kj, 5j < 0.02. For Aj = 0, the self-dual and therefore critical case, we 
have 



n sin(7rj(2m-l)/iV) 
KjC = Si 2Nsin(nj/N) hg 



sin ( (m , 



i)*/N 



sat ( (m — |)7T 



(2.29) 



The curves of 1/Kj versus A& are symmetric with respect to the vertical axis. 
2.4. Boundary Conditions and Interfacial Tension at T= 



In the study of the interfacial tensions in model systems,^E^H23c^rE§) various 
boundary conditions and different orders of taking the limit are being used. 
We shall here examine the resulting differences, using the chiral clock model 
as an example. 

In almost all of the numerical studies of the chiral clock model, fixed 
boundarv conditions are preferred. Specifically, in the work of Yeomans and 
Derrida,© they consider a lattice with C rows and M. columns, as shown in 
Fig. 4, and demand that the spins on the boundary rows have fixed values: 
a(m, 0) = r and a(m, C) = 0; but they impose periodic (or cyclic) boundary 
conditions on the utmost left and right columns. They choose to have finite 
£, but M. — > oo; that is, in the direction along the interface, the system is 
infinite from the start. 

Huse et al., in their low-temperature analysis ( -cj- ) of the wetting transition 
for the symmetric case with W = W, consider a lattice oriented diagonally 
with L rows and M columns p as shown in Fig. 5. The spins in the top row 
have fixed values r and in the bottom row they are 0. Moreover, they find it 
convenient to pin the diagonal interface in the middle by further demanding 
the spins in the upper halves of the boundary columns to have the fixed values 



5 We have made a trivial reflection. They have a "vertical" interface, whereas we have 
a "horizontal" one. 
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Fig. 4. Interface at zero temperature in a lattice with L rows and M. columns. The 
interaction between a vertical nearest-neighbor pair of spins, w n and uj n , is £{n — n'), 
while it is £(n — n') for a horizontal pair. 

r, and in the lower halves fixed values 0. That is a(0,£) = a(M,£) = r for 
< £ < \L and a(0, £) = a{M, £) = ttioi \L < £ < L. Their analysis is done 
by taking L — > oo first, then M — > oo; that is, in the direction perpendicular 
to the interface, the system is infinite from the start. 

In the analytical works of Baxter,^ ~E^ such a diagonally oriented 
lattice is being used for computational convenience, because the diagonal 
transfer matrices form commuting families. Just as in the Ising model, where 
Onsager obtained the interfacial tension© by imposing antiperiodic and 
periodic (N = 2) boundary conditions, here skew boundary conditions are 
imposed on the top and bottom boundary "rows," that is a m ,L+i = °m,o — r > 
(r— 0, ..., N — 1), for the two boundary spins in the same "column," while 
cyclic boundary conditions are imposed on the two utmost left and right 
boundary columns with a(0,£) = a(M,£). Due to such skew boundary 
conditions, a "horizontal" (actually, diagonal) interface occurs. The skew 
boundary conditions do not affect the commutation properties.^ 
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Fig. 5. Interface at zero temperature for the symmetric case with W = W in a lattice 
with M diagonal columns and L diagonal rows, with skew boundary conditions imposed 
on top and bottom rows and periodic boundary conditions on the utmost left and right 
columns. 

2.4.1. Interfacial Tension at T= 

At zero temperature only ground states contribute to the partition function 
Z, and to its logarithm, which is proportional to the free energy F. We 
would typically expect a low-temperature behavior 



- log Z = 



E B 



knT knT 



£o + o(T), 



(2.30) 



where E g is the ground-state energy, S is an entropic term related to the 
ground-state degeneracy, and o(T) stands for (usually exponentially) small 
corrections. In this subsection, we shall concentrate on the ground-state 
energy and surface or interface corrections to the bulk behavior of it. 

In the ferromagnetic case, £ (r) < £(0), £{r) < £(0) for r ^ 0, and with 
the fixed boundary conditions mentioned above, the ground state has to have 
a seam or interface. For the lattice shown in Fig. 4, the excess free energy 
at T = 0, or the increase in the ground-state energy due to the mismatch 
of spins in two adjacent rows, divided by the interface length Ai, is defined 
to be the horizontal interfacial tension e r . Hence the interfacial tensions at 
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T = are given byH 

e r = 58{-r) = £{-r) - £{0) > 0, for r = 1, . . . , iV - 1. (2.31) 

Now, from (2.27), we find that the N — 1 interfacial tensions of the symmetric 



integrable chiral Potts model at T = satisfy 

7T 

e r = 2r sin — = rt\ = e r _i + ei, (2.32) 

which is a consequence of the special form of the W weights. 

Hence, at zero temperature, this not only signifies a wetting transition, 
but for iV > 3 a more special phenomenon is taking place which we shall call 
"superwetting," with the maximal amount of interface degeneracy as each 
interface of type r is free to break up into two interfaces of types j and r—j, 
for any j between 1 and r — 1.0 It is as if we have energy levels given by spin 
operators S z for 25 + 1 = N, similar to what happens in the superintegrable 
chiral Potts model. (E§E|) 

It should be obvious that the horizontal couplings £ have no role to play at 
zero temperature. Particularly, for the three-state chiral clock model, there 
is no difference in the horizontal interfacial tension between the symmetric 
case with £{n) = £{n) and the Ostlund-Huse asymmetric case with the same 
Ai ^ but A"i = 0. 

On the other hand, if we interchange £{n) and £(n), then the incremental 
free energy due to the mismatch is now, for the Ostlund-Huse N = 3 case, 

S£(-l) = 8£{1) = 3, S£(r) = £{r) - £(0), (2.33) 

different from 

^( T l) = 3cos(yA 1 )+ v / 3sin(yA 1 ) = 2 v / 3sin(^(l + 2Ax)). (2.34) 



6 Several statements in this subsection also apply to nonzero temperatures, provided we 
replace excess energies by excess free energies per "surface area" (length) d. The resulting 
interfacial tensions or surface tensions are independent of the choice of the ensemble and 
its corresponding thermodynamic potential, as long as, with the "surface volume" (area) 
V s , each surface order parameter or its corresponding surface field vanishes.® 

7 Interfacial wetting can occur only in systems with three or more bulk phases (N > 2), 
contrary to surface wetting which can occur in the Ising model (N > 2).*EZP However, 
"surface superwetting" also requires N > 2 and interfacial superwetting N > 3. 
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Thus the excess in free energies at T = calculated for the asymmetric case 
are different in the two directions. This means that the interfacial tensions are 
anisotropic. We denote the angle- dependence by letting e r = e r (</?), whenever 
confusion may occur, with e r (0) denoting the horizontal, e r (^7i) the vertical, 
and e r (^7i) the diagonal interfacial tension. 

The diagonal interfacial tension at zero temperature can also be calculated 
by considering the incremental ground-state energy due to the mismatch of 
bonds in a lattice, for example as shown by dashed lines in Fig. 5. Following 
Baxter,^© we calculate the incremental energy per horizontal and vertical 
bond pair. For the symmetric case with W = W, we find 

e r = 25£(-r) = 4r sin(7r/iV), (2.35) 



which is double the amount in (|2.32|) and the same whether cyclic or free 



boundary conditions are imposed on the left and right boundary columns. 
Again this shows that the integrable model is at a superwetting transition at 
zero temperature. 

For the asymmetric cases, however, the diagonal interfacial tensions at 
T = are given by 

e r (lir) = 5£(-r) + 5£(-r). (2.36) 

For the three-state Ostlund-Huse model at Ai = |, we have 

SS(-l) = V3, 8E{-2) = 2a/3, SS(-1) = 6E(-2) = 3. (2.37) 



Substituting these into ( |2.36|) , we find that |ei(|7r) is greater than ei(0) 



given in ( |2.32j ) and ^62(^71") is smaller than €2 (|vr) in ( |2.32j ). Hence this 
diagonal interfacial tension, as well as the vertical interfacial tensions, in the 
Ostlund-Huse asymmetric case with W 7^ W, do not satisfy the superwetting 
condition even at T = 0. 

At zero temperature, we can also easily find the interfacial tensions for 
general angle < <p < ~ir, as 

e r ((p) = N{(p)- 1 (8E{ r r) cos^ + 8£{-r) sin^), (2.38) 

which can be normalizeda per horizontal (or vertical) bond using 

Af(<p) = maxfcos(y9),sin(y9) J. (2.39) 

8 For the excess energy per unit length we would have to use N{(p) = 1 . 
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Such a ground state is highly degenerate as there are many configurations 
of the interface with a given number of horizontal and vertical bonds. For 
ground-state wetting we need that either a horizontal or a vertical piece 
of the interface wets. So, for N = 3 the system is not wet when both 
< Ai, Ai < -r, and it is wet when at least one of the Ai, Ai > ■? within the 
interval < A]_, Ai < | . A ground-state wetting transition occurs for 

A x = I, < Ai < \ or < A x < 1, Ax = ±, (2.40) 

which is the boundary of these two regimes. 

At Ax = f, we find from ( ggg ) and (|2T34l) that 



<S£(-1) = 0, <J£(-2) = 3, <J£(-1) = <J£(-2) = 3. (2.41) 

Now we can use ( |2.36| ) to find that the diagonal interfacial tensions do satisfy 
the condition for the onset of wetting: 62 = 2e\. This means that the wetting 
transition of the diagonal interface of the Ostlund-Huse model occurs at 
T = and Ai = |. It is interesting to note that the chiral melting line starts 
at the same point. The finite-strip calculations of Yeomans and Derrida® 
show that even the vertical interface (which is parallel to the chiral field) is 
wet at this point. It is interesting to investigate whether the wetting curve 
of the diagonal interface is identical to the chiral melting curve. 

We shall now consider some other subtleties of the interfacial tensions, 
as they relate to different boundary conditions and different orders of taking 
the thermodynamic limit. We shall restrict ourselves to diagonally oriented 
lattices only. 

2.4.2. The Limit M^oo, then L — ► 00 

If one lets the number of columns M — ► 00 first and the number of rows 
L — > 00 afterwards, then the boundary condition imposed on the columns 
should not have any effect. This can also be seen by considering the elements 
of the column transfer matrix T^ with finite L -- dividing by the largest 
eigenvalue of T^ 00 - 1 (the bulk term). At non-zero temperature, these elements 
are all positive and thus the Perron- Frobenius theorem holds. Consequently, 
the largest eigenvalue is nondegenerate. In the limit M — > 00, only the 
largest eigenvalue survives. This shows that the resulting interfacial tension 
is independent of the boundary conditions imposed on the columns for T > 0. 
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Nevertheless, in this limit, the interfacial tensions depend heavily on the 
boundary conditions imposed on the top and bottom rows. We may consider 
the interface (or domain wall) as a random walker,® who tends to walk in 
the direction of higher probability. For the skewed boundary condition, this 
allows an interface winding around the cylinder of length L and perimeter 
M crossing the seam of modified bonds several times. Thus for the Ostlund- 
Huse model, with© 

x = w x = W(-1)/W(0), y = w 2 = W(-2)/W(0), 
z = Wi = W(-1)/W(0) = w 2 = W(-2)/W(0), 

x < z < y, (2.42) 

the ei interface prefers to walk perpendicular to the chiral field and then 
continue by crossing the seam, while the e 2 interface tends to walk parallel 
to the chiral field, as shown in Figs. 6a and 6b. Thus one does obtain the 
horizontal interfacial tension ei(0) and the vertical interfacial tension e 2 {^) 
using a diagonally oriented lattice with skewed boundary conditions. For 
fixed boundary conditions on the boundary rows, one obtains instead the 
diagonal interfacial tensions e r (j7r). Specifically, we find 

skewed b.c: lim lim k^T \og{Z r / Zq) = min e r (8), 

fixed b.c: lim lim kvT \og(Z 0r / Z 00 ) = e r (j7r), (2.43) 

where Z T and Z$ r denote the corresponding partition functions. 

In the latter case, with fixed boundary conditions for the boundary spins, 
we have free open boundary conditions for the domain walls, which can touch 
but not cross the boundaries. On the other hand, a seam due to skewed 
boundary conditions can be moved to any place using gauge transformations. 
In spite of the Perron- Frobenius theorem, these two cases are very different, 
even in the limit M — > oo, when the size of the row transfer matrix becomes 
infinite. 
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(b) 




Fig. 6. (a) For the case with skewed boundary conditions in both directions, an interface 
is allowed to wind around the torus of vertical perimeter L and horizontal perimeter M , 
crossing the seams of modified bonds several times. Thus for the Ostlund-Huse model, 
a single-step interface with tension t\ prefers to walk perpendicular to the chiral field 
and then continue by crossing the seams, (b) A double-step interface with tension 62 tends 
to walk parallel to the chiral field, (c) If fixed boundary conditions are imposed on the 
top and bottom rows, and free boundary conditions imposed on the outer columns, it is 
energetically more favorable for the system to arrange itself into a configuration with a 
mismatch. Therefore if we let M — ► 00, and fix the ratio M/L > 1, the excess in free 
energy at T — can have any arbitrary value depending on M/L. 



2.4.3. The Limit L^oo, then M^oo 

If we take the limit that the number of rows L — > oo first, then it is necessary 
to pin the interface roughly at the middlea in order to get a true e r interface. 
Otherwise, because there are L positions to place the e r interface, while 
there are ( L ) possible positions to place the r e\ interfaces, the term with r 
t\ interfaces dominates the partition function, in this limit. 

Now it is easy to see that the boundary conditions on the top and bottom 
rows can have no impact, as an interface of finite length M could not possibly 
reach the boundary rows which are infinitely far away. On the other hand, 
with one end of the interface being pinned to the middle, cyclic boundary 
conditions imposed on the boundary columns force the interface to come 
back giving the diagonal interfacial tension, while free boundary conditions 
on the boundary columns allow the interface to wander and settle to its lowest 
energy configuration, yielding the minimum of the interfacial tensions. 

2.4.4. The Limit L,M^oo, with fixed M/L 

In exact calculations, however, it is cumbersome to pin the interface in the 
middle. As Baxter® chooses to let M, L — > oo simultaneously, we shall 
examine the case when they are proportional to one another. Since the 
Perron-Frobenius theorem may not hold on infinite matrices, a difference in 
boundary conditions could play an important role. 

The full complexity seems to arise when M/L > 1. It is most interesting 
to consider the case when fixed boundary conditions are imposed on the top 
and bottom rows and free boundary conditions on the columns. Then it is 
easily seen that it is energetically more favorable for the system to arrange 
itself into a configuration with a mismatch as shown in Fig. 6c. Therefore, if 
we let M — > oo and fix the ratio M/L > 1, the excess in free energy at T = 
can have any arbitrary value depending on M/L. A similar problem arises, 
when skewed boundary conditions are imposed on the top and bottom rows. 
For if the ratio M/L is an integer, or free boundary conditions are imposed 
on the boundary columns, we will obtain the minimum of the two interfacial 
tensions shown in Figs. 6a and 6b; otherwise, the interfacial tension will be 
a function of M/L. 



9 This is equivalent to the subtraction procedure used by Baxter,® who needs to omit 
an L 2 Z\ term from Z 2 in order to obtain e 2 , see his (5.27), (A9), and surrounding text. 
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For M/L < 1, however, this complexity disappears, and cyclic boundary 
conditions on the columns give the diagonal interfacial tensions, whereas free 
boundary conditions give the minimum of the interfacial tensions over all 
directions. Thus by taking the limits in an appropriate way, we may calculate 
horizontal or vertical interfacial tensions using a diagonally oriented lattice. 

2.4.5. Single Interface at Low Temperatures 

Let us now discuss the effect of boundary conditions on the restricted 
partition functions with one domain wall. As before, M denotes the size 
of the system in the direction of the interface, and L the size perpendicular 
to it. Besides the partition function without the domain wall Z (L,M), we 
can introduce various partition functions with a single domain wall. First, 
we can introduce Z lpin (L, M\s) for the partition function with the interface 
pinned in the middle on one side and pinned s steps from the middle on the 
other side. We can also introduce Z lpcT (L,M) for the periodic case, where 
the interface starts and ends at the same, but otherwise free, position, and 
Z 1 f ree (L, M) for the case where the interface is left free on both sides. These 
three partition functions are closely related when L becomes large, namely 

Z 1|PCT (L,M) Z 1|Pin (L',M|0) 

Z (L,M) ~ l™oo Z (L',M) ' 

Zi |fr ee(£,M) +2° Z 1|Pin (L',M| a ) 

Z (L,M) ~ ^i'-oo Z (L',M) ' K ' J 

Therefore, these quantities are easily evaluated at low temperatures, for 
which overhangs can be ignored, using e.g. random walks, transfer matrix 
techniques, or Szego's theorems for Toeplitz matrices. 

For the horizontal interface as in Fig. 4 and using ( |2.42| ), we find 

1 Z 1Mc (C,M) (l-xy)x 

hm hm — -log— - 1 — — — — — ~log- 2.45 

M^oc/:->ocM £Z (£,M) & (l-x)(l-y) v ; 

and 
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1 Z 1>pcr {C,M) 

lim lim — - log — — - — — — 

M^ooC^ooM CZ {C,M) 

1 Zi |flx (£,A4) (l-xy)x 

w lim lim -r-7l°gT7 J — ' , ' & log- —, (2.46) 

£ _ooA<-.ooA^ 5 Z , fix (£,.M) & (l- v ^y) 2 ' V ' 



where "fix" stands for various boundary conditions which keep the interface 
continuing within a horizontal strip of width C The results differ only by 
exponentially small terms in the temperature T, following a leading term of 
order 1/T, showing that "kinks" are important. 

However, for the diagonal interface as in Fig. 5, we find 

lim lim -*- log ^j^^ « 2 log(x + x) (2.47) 



M-ooL-kxjM LZ (L,M) 



and 



lim lim — log — L ^— - — - — - 
M-»oo£-»ooM b LZ (L,M) 

« lim lim i log I^T 1 ^ w 2 lo g( 2 ^)> (2-48) 

where the factor 2 is to have agreement with Baxter's convention. Now the 
differences occur in the order 1, and the corrections are entropic in nature. 
From the results ( 2.45 ) to (|2.48|) , we see that the effects of the boundary 



disappear when the system is reflection symmetric with respect to an axis 
perpendicular to the interface, so that x = y in (|2.45|) and Q2.46|) , or x = x in 



( |2.47| ) and ( |2.48| ). In those two cases, for which the interface is perpendicular 



to the chiral direction, we can use ref. 49 to obtain further detail. 



2.5. Wetting Transitions 

Wetting transitions can occur when there are different types of domains 
and domain walls. More precisely, the wetting transition is defined as "the 
bubbles of B domain absorbed on the A\C interface merge into an essentially 
macroscopic layer of B domain which wets the entire interface," (see e.g. page 
377 of ref. fJ5p . We need to compare a configuration with a domain wall A\C 
with a configuration with two domain walls A\B and B\C. Hence, the wetting 
temperature T w is defined by 

T w = min T w (0), with T w (0) from e 2 (0, T) = 2 ei(0, T). (2.49) 

6 
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As the interfaces can be oriented differently, the wetting transition can 
occur at different temperatures T w (9), but its minimum is the true wetting 
temperature for given chirality A x . 

We note that we need to compare the interfacial tensions 62 (#) and 
2ei(#) for the same angle 6, as is implicit in all the calculations mentioned 
above. (c^rEiP It is often energetically more favorable to wet than to turn the 
interface through an angle for which it is necessary to flip macroscopically 
many spins. We should be careful not to calculate T w from mine2(0, T) = 
2 minei(#,T), as a calculation with skewed boundary conditions in both 
directions could lead to and which would cause us to overestimate A wet (T). 

It seems that an interface perpendicular to the "chiral field" Ai wets 
first. In fact, Huse et al. calculate the wetting curve by considering diagonal 
interfaces for the symmetric case, with A ver = Ah or , and horizontal interfaces 
for the Ostlund-Huse model with A ver 7^ and A hor = 0. Particularly, for 
the diagonal interface of the Ostlund-Huse model, its wetting line would 
have to start at zero temperature at Ai = |, above the wetting curve of 
the horizontal interface, which wets at Ai = | and T = 0. We note that 
a diagonal interface is a superposition of many allowed walks with a given 
number of horizontal and vertical bonds. So with a chiral field in the vertical 
direction, horizontal segments of the interface will wet at 7, while vertical 
segments will not wet, lowering A wet (T = 0) from | to \. 

It is easy to extend the low-temperature analysis of Huse et al.tilf to the 
case with K\ 7^ K x and one finds that, at low temperature, the wetting line 
for the Ostlund-Huse model, (x = y = z) , is given byE3 



y 



X 



6K1 sin (±vr(l - 4Ai))] = 1 + 2z 2 + 0(z 3 ), z = exp(-6ifi), 

(2.50) 
while the integrable line is given by 

^ = l + 2z + 0(z 2 ). (2.51) 

.X. 

The integrable line, denoted by line 3, and the wetting line, denoted by line 4, 
are shown in Fig. 7b. They rise faster than any power due to the importance 
of kinks, see also the text below 



10 There is a small discrepancy with a factor 3 given by Huse et al.© instead of 2 as 
given below, which we believe to be due to a misprint. 
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Fig. 7. (a) Intcgrable lines in the 3-state chiral clock model. The solid curve is for the 
Ostlund-Husc fully asymmetric case and the dashed line is for the symmetric case. The 
crosses are the numerical results of Stella et al. for the fully asymmetric case, (b) We 
enlarge the part of the graph for T small and Ai near i, with the intcgrable lines marked 
as lines 1 and 3. We have also plotted the results of Huse et al. with line 2 denoting the 
wetting line of the symmetric lattice, and line 4 the one for the fully asymmetric case. 
The graphs show that the intcgrable lines are near but below the wetting lines, with the 
largest relative deviations in the low-temperature regime. 
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For the symmetric lattice with W = W, (x = x,y = y), however, we find 
that the diagonal interface wets first, with the wetting line given by Huse et 
al. as 

(2.52) 



L = 3 4 

2 ° ' 



x* 2 
whereas the integrable line is given by 



y 



2 + -- 



(2.53) 



These are also plotted in Fig. 7b as lines 2 and 1 respectively. The wetting 
curve for the horizontal or vertical interface in the symmetric case also starts 
at Ai = | and T = 0, but rises in a much steeper way. Since all the interfaces 
wet on this curve, it would be interesting to compare the behaviors near this 
curve and near the chiral melting curve of the Ostlund-Huse model. 

We show in Fig. 7a the integrable line for the 3-state chiral clockt^rEi' 
model. The solid curve is for the Ostlund-Huse asymmetric case, while the 
dashed line is for the symmetric case, with equal horizontal and vertical 
interactions. We can see that both curves end up at the same point for zero 
temperature. The numerical results of Stella et a\. < & for the wetting line in 
the Ostlund-Huse model are represented here by crosses. 

In Fig. 7b, we enlarge the graphs for T small and A x near j, with the 
integrable lines marked as lines 1 and 3 added to the wetting line results 
of Huse et al. with line 2 for the symmetric lattice and line 4 for the fully 
asymmetric case. The curves show that the integrable lines are near but 
below the wetting lines. For < T < T c , therefore, contrary to what we have 
conjectured before,^ and in agreement with Baxter,®© we find that the 
integrable model is in the non-wetted region. 

2.6. Low- Temperature Regimes 

As suggested to us by Baxter Jill we can approximate the Boltzmann weights 
of the integrable chiral Potts model, 



w„. = 



W{-n) 
W(0) 



n-l 



7V 3=0 



sin(jA - 4>) 
sin(jA — 6) 



11 Private communications. 
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w n = 



W(-n) 



n— 1 

n 



W(0) A" f = L [sm(jX-e) 



sin(jA - 0) 



where 



A- n 



+ 



0-g \ N - sin W) -^N^sm{N<f)) 



sm(N6) 



sm(N9) 



(2.54) 



(2.55) 



see also ( p.7|) , (|2.21|) , and section [|, by replacing 6 1 , #, 0, inside the products 
in ( 2.54J ) by their zero-temperature limiting values. For the symmetric case, 
with 

-\<9 = 9<-l\, = = # + ±A, (2.56) 

we can express 9 and in terms of k using (|2.25|). We find successively 



k = cos(2N9), sin(N9) = -\j~{l - k), sin(iV0) 

1 + ife , -r 



|(1 + k), 



A = A 



9 = 9 



k 
~k 



u/ ' 



1 



._ ^ _ arcsin ^1(1 _ k)) = - — . 

The zero-temperature limit then corresponds to 



— arcsin ^±(1 + £;), 

TV 



9 = 9 



-A, 






oo, 



and we find the asymptotic low-temperature formula for the weights 

_((i-|M) 

3=1 



W, 



U.K 



i n sm 
A" Li sin(jA) 



(2.57) 



(2.58) 



(2.59) 



Here we can allow A ^ A, as this corresponds to a gauge transformation. 

Baxter applied a Bethe Ansatz method with skew-periodic boundary 
conditions,® shown in Fig. 5, together with a complicated subtraction 
method, in order to obtain the first two diagonal interfacial tensions in the 
low-temperature regime. Indeed, substituting 



/ = 9 = tan UX), Xo = 0, 



(2.60) 



corresponding to the symmetric case, his results (5.36) and (A12) reduce to 
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fl 



k B T 


- -iug 


C2 

k B T 


2 £l 
k B T 


or, equivalently, 






ei 




fc B T 




62 



- log^Wi) = 2 log (A cos(|A)) , 



log((l-/ 2 ) 2 ) =21og(l-tan 4 (iA)), (2.61) 



2 log A + log(l + tan 2 (±A)^ 

41ogA + 21og(l-tan 2 (±A)). (2.62) 

These results are useful to verify that the exact calculations© based on 
analytic continuation from the superintegrable case are correct. 

For the symmetric case, there is an easier Bethe Ansatz calculation than 
the one done by Baxter,® namely one with free boundary conditions. This 
gives the desired results directly, without a subtraction procedure. As we 
do not only have to consider "collisions" of interfaces, but also reflections 
at the boundary, we use an extension of the Bethe Ansatz^ 2 ! similar to 
the ones introduced by Gaudin.® But as this requires further consistency 
conditions, our method does not apply to asymmetric cases that differ from 
( |2.59|) by more than changing A to A for w. We shall outline our Bethe Ansatz 
calculations in section ^. These calculations show some special features on 
the integrable line, which may be related to some pre-wetting phenomena. 

We note that the Bethe Ansatz method to calculate interfacial tensions 
in the low-temperature limit is similar to the Muller-Hartmann-Zittartz 
approximation® in ignoring overhangs. Moreover, that method can be also 
used for non-symmetric and non-integrable cases, with the condition 

y / x 2 = e 6^ 1S in( 7 r(l-4A 1 )/6) ) y/ -2 = ^K, sin^l-^/e) finite _ (353) 

However, for more than two interfaces the additional condition y/x 2 = y/x 2 
is needed for the Bethe Ansatz to work. 

The above procedure ( p.58|) to obtain the low-temperature expansion also 



works for asymmetric cases with 

0,0->-A, fc->l, 0=-A-0, with 074-AorO. (2.64) 

Later on in section f| we shall see that this represents just one part of 
the integrable manifold in the phase diagram for the 3-state chiral clock 
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model. So we shall call models with ( |2.64| ) "quasi-symmetric" cases. It is a 



regime which does not include, nor even border, fully asymmetric cases as 
the Ostlund-Huse model, for which the vertical weights satisfy the condition 
W(n) = W(—n) while the horizontal weights are chiral, W(n) 7^ W(—n). 

As we shall show for N = 3 later in section f|, to have W(n) = W(—n), 
we must require 

= _A-0, with A = l. (2.65) 

Together with the integrability condition (|2.55|), we find then that there can 



be only two free variables, say 9 and 0, such that 

e = -\-\{e-<f>i f=\(e-<i>), (2.66) 

Because of that, in the low-temperature limit given by 9, 9 — > —A, we must 
also have — > — A. In order to see whether the integrable model is or is not 
in the wetted region, it is necessary to include higher order terms in the low- 
temperature expansion. In fact, for the three-state Ostlund-Huse model, the 
condition K = K reduces the free variables to just one, and in the low- 
temperature limit, — > — |7r + 50, we find 9 = — |7r + 50 p with p = (1 + v3). 
(The details are delegated to section f|). Therefore, we have 

— = 0(50) ~ w 1 = w 2 = z = , n 1 -rr. (2.67) 



Wi z x z sm[|7r — ^(0 — v)\ 

Because these are of the same order of magnitude, we have to take kinks in 
the interface into account,^ as also noted below ( |2.46| ). 

It seems that the largest asymmetric integrable regime is described by 
— > 0, 9, 9, — > A. For N = 3, this corresponds to K = K and A = pA, 
{pj^O, 1), which again reduces the parameters to just one. We find 

= --7T + 50, 9 = --7T + 59, = -50, 9 = --7T + 50, (2.68) 

with 

/ 2 \ Vb 2 +b+i _ _ i \ 6-i 

50^50 I ^50 , 50^50, 50 « 50 I -^50 I , (2.69) 



where 



b = sin (1(1 + p)ir) /sin (1(1 - p)ir) . (2.70) 
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For limits as in ( |2.64| ), we shall show in section |] that 

Ik it- 1 k , 

A 1 = + -.., Ai = --^ + ---, 2.71 



5. ±ms snuws mat uuin z_v aiiu 

as T — >■ 0, which is rather restricted. 



with k and k some positive constants. This shows that both A and A . 



2.7. Baxter's Exact Interfacial Tension Results, Symmetric Case 

In two recent papers,®® Baxter obtained several new exact results for the 
interfacial tensions in the integrable chiral Potts model. At this moment, it 
is not clear to us if or how we can extract from his results the interfacial 
tensions for the fully asymmetric Ostlund-Huse case. Here we shall work out 
his results in explicit detail only for the fully symmetric case, leaving the 
results for the more general case implicit for now. 

In equation (42) of Baxter's second paper, we can substitute® 



«* = *« = ^r 


IK A q = ] J i _ k , Vp = l, 


Vq =V, 


A 1 


ll h' 




N/2 Iv q l 


/ 1 + *. ^VT^, m - 


= m = 1 



(2.72) 



resulting in the very explicit formula for the interfacial tensions 



€r o ( \ 8 Pa sin (f) f , h N -y N , 0T * 

— — = 2v r m = - / dy—— f— — — artanhW ^^, 2.73 

k B T n Jo 1 + 2y cos(^) + y z V 1 ~ r] n y n 

in the symmetric case. Here rj is a temperature-like variable on the integrable 
curve, defined in (|2.72f ). 



In the low-temperature region, 

k' « |JV(1 - rj) -> 0, (2.74) 

we can expand (|2.73|) as 

e 2r 

~ logd^-a^ + O^'), (2.75) 



k B T N 
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where the constant is given by a dilogarithm integral, which seems different 
from the ones studied recently by Kirillov and many others,^® 



4 yl 

Cr,N = — / 
TV JO 


dx Sln(f) 

l + 2xcos(f)+x 2 


Hl-z" 


Vcos 2 (|rA) 


n 


cos 4 Q(r 
cos 4 (§( 


-2j + 1)A)n 

r-2j)A) / 


This also satisfies the sum 


rule 










JV-l 

E* 


3,N z 


= log JV. 





(2.76) 



(2.77) 



We have not found a direct way yet to prove the last equality in (|2.76|) , 
except for the cases N = 2 (Ising) and N = 3. In section |6], we shall outline 
an indirect proof using a Bethe Ansatz method. We originally guessed the 
formula by performing numerical integration to 100 places and expanding 
the exponentials of the integrals in periodic continued fractions, which gave 
us 

_ _ 4 _ 9 

o-i,2 = log 2, 0-1,3 = log-, ^2,3 = log-, (2.78) 

o M = log(2(2 - \/2)), o 2 , 4 = -log(2(2 - yftf), 

3 . 4 = log(4(2-V2)), (2.79) 
o li5 = log(2(5 - V5)/5), o 2>5 = log(5/4), 

03.5 = log(32/(5 - Vbf), o 4i5 = log(5(5 - Vbf/10). (2.80) 

Combining these results with Baxter's two Bethe Ansatz results© for r = 1 
and 2, see also (|2.62j ) in the previous subsection, we originally guessed Q2.76Q . 



For the critical region, 



iiV ^ 1 



v ^ « ±k -> 0, (2.81) 



we can expand (2.73) by first expanding the artanh in a Taylor series of its 



argument. Then, it is straightforward to rewrite the integrand in ( |2.73| ) as a 
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power series of y and rj, multiplied with Ji] N — y N . Each term of the series 
leads to a Beta function integral, B(x,y) = V(x)T(y)/T(x + y). We find 

e r _ 8 sin(f )B(±, \) i N+1 _ 8 sm(^f)B(^, |) i N+2 + Q 1jv+3\ 



A; B T vr(iV + 2) ' vr(iV + 4) 

rEl (2 ' 82) 

giving both the leading term© and the first correction term. Coefficients of 

further terms r]2 N+ \ for j = 3, 4, . . ., can also be obtained easily. 

2.8. Other Exact Results and Further Speculations 

A great deal of progress on the chiral Potts models has been achieved. The 
exponent a for the specific heat — not to be confused with the independent 
variables for the integrable weights given in (|2.13|) - - has been obtained 
by Baxter^ for the case with real positive Boltzmann weights and the 
exponents (3j of the one-point functions < cr J > have been conjectured by 
Albertini et al.® They are 

«-l-^, Pi- m2 ■ (2-83) 

In fact, Albertini et al.® also conjectured the exact low-temperature formula 

< a J >= (1 - k l2 f = k 2 ^ , (2.84) 

which generalizes the Onsager result^Ho for the spontaneous magnetization 
of the Ising model. 

For N = 3, we find a — |, which is identical to the three-state Potts 
model result. It is generally believed that the chiral field Ai is a relevant 
variable, and the free energy in the scaling region can be written as 

F(T, A b A x , H) = F(T, A 1 ,pA 1 , H) = \t\ 2 ~ a X(gA 1 /\t\t hH/\t\% (2.85) 

for some crossover exponents 4> and ip and with t = T/T c — 1, where T c is the 
Potts critical point. The ratio of the two chiral fields Ai/Ai is believed to 
be an irrelevant variable, which may only change the two constants g and h. 
The exact calculation of BaxterW gives 

F^A^A^O)-^ 2 ! 2 "", (2.86) 
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where k is the elliptic modulus given by ( |2.23| ), with k = for T = T c and 
Ai = 5i = 0. 

For k ~ 0, we find the integrable line lies on the curve given by 



-9.959616982 (1 + p 2 )A\ + 0(t 2 , tA 2 ± , Af), p= Ai/Ai, 



and 



108 (2 + V3) Kl 



■t = -4.981050242 L 



(2.87) 



(2.88) 



1 + K lc (7^3 + 12) 

with Ki c = Ji/k&T c = |log(l + \/3), see also eq. (|4.26|) . In the above 
equation, p — 1 corresponds to the symmetric case and p = to the Ostlund- 
Huse asymmetric case. The analytic expansion Q2.87 ) for the integrable line 
is similar in form to the usual expression for the nonlinear thermal scaling 
field© 

i=t + c 2 A\ + • • • , (2.89) 

in terms of which the free energy can be rewritten as 



F(T,A U A U H) = \tr a X( j gA 1 /\t\*,hH/\t 



(2.90) 



But to say more about this requires more detail from exact results than we 
currently have available. 



From the most recent results in (|2.82|) , we conclude that we have the 
following asymptotic expansions for the surface tensions on the integrable 
line, 



e r = \t\? + *D r (\t\*), 
which we must compare with the forms required by scaling 

e r = It^A^Ai/ltl*). 



(2.91) 



(2.92) 



Now, as the temperature is an analytic function of Ai on the integrable line, 
we need 

Ai ~ |t|3, (2.93) 

and we find 

(i=~ + —, <j> = ~-±=. (2.94) 



1 1 1 

The critical exponents are most conveniently expressed in terms of the 
scaling (conformal) dimensions Xj and jjj, indicating how a local density (or 
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order parameter) rrij and its corresponding field scale with a typical length 
scale R in the problem, 



m 



■j ~ R~ x i, hj ~ R~ y i, Xj + yj = 2. (2.95) 



In our model, the subscript j takes the values T (thermal), 1, . . . , N — 1, A 
(chiral), with hi = t = (T — T c )/T c , rrij =< a^ > for j = 1, . . . , N—l, and 
h A = Ai. 

In terms of these, the usual critical exponents are 



1 

V = \k= — , 


Tjj ^%ji 


2-a= — , 


^ = ~ 


V? 




Vt 


Vt 


Pi = % 
Vt 


_Ui- Xj _ 

ho 

Vt 


£j Jb •£ Jb j 

Vt 


Xj 



where we use ipj for the gap exponents (or crossover exponents), as the more 
conventional notation A would lead to confusion. We are writing = ip& for 
the chiral crossover exponent in ( |2.85|) . 



This leads to the scaling dimensions 

4 j(N-j) N + 6 

XT ~ iV"+2' Xj ~ N(N + 2)' XA ~ iV"+2' 

2N N-2 . „, 

Vt = ]yq^» Vi = 2 ~ x i' V* = J^^ =yr-l. (2.97) 

For N = 3, we find = |, xa = I, in agreement with earlier predictions. ^& 
For general N, y& = Ur — 1 provides the chiral exponents of the Fateev- 
Zamolodchikov model, ca* 

The most detailed results have been obtained for the superintegrable 
model with weights given by 

W(n) ((l,q)) 0|W W(n) ((l,/?/a;))o,n 

W(0) " ((1, /?))„,»' W(0) " ((1, <*))„,„ ' { ' } 

For the corresponding hermitian quantum spin chain, Albertini et al.feP 
have calculated results for the excitation spectrum, demonstrating the special 
role of level crossings, as the Perron- Frobenius theorem does not apply to this 
case. Baxter® has calculated the bulk and surface free energies, horizontal 
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and vertical interfacial tensions and finite size corrections for the lattice 
shown in Fig. 6 with free or fixed spin configurations on the left and right 
columns. He has found that the critical exponents for these physical 
quantities are given by 



« = l-f, 


a =2-2- 

«s ^ N , 

_ 2 _ _ 

-~ at) /* ver ^ver 


l^hor ^hor ' 




^hor + *Arer = 2 - <X 



(2.99) 

These are extremely interesting and puzzling results, deserving a great 
deal of attention. Even though the weights are asymmetric, one should not 
associate z/ ver to be in the direction of the vertical weight, because the lattice 
is oriented diagonally. The two different correlation lengths are obtained from 
the decay of the finite size effects of such diagonally oriented lattices. Since 
the diagonal lattice is symmetric under 90° rotation, one is forced to conclude 
the different boundary conditions give rise to the two different exponents for 
correlation lengths. 

Yet it is remarkable that the bulk properties should be influenced by the 
boundary, particularly as the partition functions in the superintegrable case 
are real. In the most recent work of Baxter JEao he imposes cyclic and skew 
boundary conditions on the two directions, finding that 2/x = 2 — a = 1+2/N. 
Assuming jjl = v, then the usual scaling relation again holds and, with 
such boundary conditions, we would have that the vertical and horizontal 
correlation lengths have the same exponents with u VCT = z/h or = 1/2 + 1/N. 

In the superintegrable case, the weights are complex. By a gauge 
transform, we can make the weights real, but we cannot make them all 
positive. Therefore the Perron-Frobenius theorem does not hold. This may 
be the only plausible reason for the above strange behavior. Because the bulk 
properties can be shown to be boundary dependent for non-positive weights 

- complex weights are typical within the integrable chiral Potts model - 
the combination of using Z-invariance properties and analytic continuation, 
as used by Baxter, may or may not lead to valid results for the positive-real 
Boltzmann weight case. For this reason, we find it necessary to verify some 
of the calculations of Baxter by low-temperature expansions. 

For real and positive weights, such boundary-dependent behavior is not 
conceivable. Furthermore, because the anisotropy p = Ai/Ai is an irrelevant 
variable, the exponent v cannot be a function of p either, therefore we must 
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have ^or = ^ver = ^dia = v. Baxter*!^!® found that \x$ xa , = 1/2 + l/N. 
Hence, assuming the scaling relation n = v, one finds that the scaling relation 
2v = 2-a = l + 2/N again holds. 

2.9. A New Model 

In |5.2| we shall present plots, for iV = 7, of l/K 2 and 1/K 3 versus 1/Ki, 



and of A2 and A3 versus Ai. The curves look like straight lines. Plots for 
N = 5, 6 also give the same impression. Yet an explicit calculation shows 
that the ratios Kj/Ki and Aj/Ai are almost but not quite constant. We 
therefore propose a new model with constant ratios, given by 

E(n) T ^ sin(7r/iV) (2it ( 



which has only two parameters. This may be good for further numerical 
studies. For A = \ and K = -E/k B T, (j^TUUp and (^27|) are identical. 
Hence at zero temperature, it goes through a superwetting transition at 
A — j. In section ^|, we shall describe some of the symmetries of this model 
and show that for A = \(N + 1)/(N - 1) we have 8(0) = £{-!). Hence, 
its ground state is highly degenerate, as it is in the chiral clock model; ^O 
thus this model, like the chiral clock model, can also be used to describe 
commensurate- incommensurate transitions. 

For T > 0, we shall show that for a certain value of A, this model is 
not very different from the integrable chiral Potts model. We can use the 
universality hypothesis to argue that this model in two dimensions, having 
the same Z^r symmetry as the integrable chiral Potts model, must have the 
same exponents. Hopefully, the exact results obtained for the chiral Potts 
model can be used to gauge the accuracy of the numerical studies. 

3. CONSISTENCY EQUATIONS AND INVERSE PROBLEM 

In this section we derive the consistency equations that the weights must 
satisfy in order to be put in product forms as in ( |2.13 ). We also express the 
variables a and j3 in terms of these weights. 
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3.1. Consistency Conditions for General JV 

Let us introduce the notation 

_ W(n)W(N-l) 
tW- W (n-l)W(0)- [6 } 

Then substituting (p. 13 ) and 02.141 ) into it, we find that the A(x) terms 
cancel out leaving 



/(«) 



(1 - au n )(l - (3) 
(I ~ f3u n )(l - a) 



1 < n < N - 1. 



(3.2) 



Now we use ( |2.20| ) to rewrite 

cot (vm/iV) 



l + u n 



UJ' 



cot 9 = —i 



l + a 



cot « 



l+P 



(3.3) 



1-a 1-/3 

Consequently, ( |3.2p become AT — 1 linear equations for cot 9 and cot 0, i.e. 

cote -f(n) cot0= (f(n) - l\cot(im/N), l<n<N-l. (3.4) 

The determinant of any three of these equations must vanish in order to have 
nontrivial solutions. This gives the consistency equations 



sin 



(n-2)7T 

N 



/(n)-l)(/(2)-/(l) 



sm 



rnr 

~N 



(/(l)-l)(/(n)-/(2) 



- (3 - 5) 

for 3 < n < N — 1. Thus if the weights W(n) and W(n) satisfy these N — 3 

consistency equations, then they can be put in the product form. 

We note that (|3.5|) is equivalent to the more general equation 

(f(ni) - f(n 2 )) (/(n 3 ) - /(n 4 )) _ sin(^(m - n 2 )) sin(-^(n 3 - n 4 ; 



(/(ni)-/(n 4 ))(/(n 2 )-/(n 3 )) sin(f( 



«-i — ri4 J sm 



■(*< 



n 2 - n 3/ 



which is solved by any expression of the form 

X\UJ n + X2 



/(n) 



x 3 ct; n + x 4 
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(3.6) 
(3.7) 



One can solve any two of the equations ( |3.5| ) , yielding 

/(n 2 )(/(ni) - l) cot(imi/N) - /K) (/(n 2 ) - l) cot(7rn 2 /iV) 



cot0 



(/(na) - /(m 



(/(m) - l) cot(7rrn/iV) - (/(n 2 ) - l) cot(vrn 2 /iV) 
cot0=^ '- -, ^ 1 3.8 

(fM - /(m) 



which expresses 9 and in terms of the weights W. 
Using the integrability condition ( |2.22| ) we find 

QQ = 1, (3.9) 

with 

IT 1 = sin(vr/iV) (cot(0 -9)- cot(n/N)) (3.10) 

and a similar equation for Q. Since we have expressed 9 and in terms of 
the ratios of weights f(n) and similarly 9 and in terms of the f(n), we can 
substitute ( |3.8| ) into ( |3.10| ); then ( |3.9| ) gives a relation relating the W and W 
weights. 

For N = 3, (|3.9|) and ( J3.10I) are the same equations that we presented 



elsewhere without derivation,M except for a slight change of notation, as K 
and A in the earlier work are changed to 2K\ and Ai here. These equations 
are useful to check whether a particular model is integrable or not. Yet even 
in the N = 3 case, we find it more convenient to use the parameters 9 and 4> 
instead. Moreover, even though we can in principle use the equations (|3.5| ) 
to calculate the Kj and Aj for j > 2 in terms of K\ and Ai and (|3.9|) and 



( |3.10|) for the relations between the Kj and Aj and the Kj and Aj, we have 



found out that these algebraic equations become more and more complex to 
solve, as iV increases. In the next subsections, we shall present an alternative 
way. 

3.2. Interaction Energy Parameters of Chiral Potts Model 

We may rewrite (|2.13|) in terms of the variables 9 and as 

W{n) /sin(iV0) \ 1/N sin(9 + nn/N) 



W{n - 1) \ sin(N9) J sin(0 + tcu/N) 
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(3.11) 



Together with fl2.2| ) and ( |2.5| ), we find 



S(n) - S(n - 1) _ W(n) 



JV-l 



fc R T 



W(n-1) 



A + £ n = X! 



fjW' 



j a 



i=i 



with 



.4 



1 sin(iV0) 
iV ° g sin(iV#) 



-B„. = lot 



sin(0 + Tin/N) 
sin(0 + nn/N) ' 



knT 



{l-u,->). 



The Fourier coefficients Tj in (|3.12 ) can now be written as 



A' 



N 



Nr = NA + Y,B n = 0, 



t j 



N- l Y.^ n3 B ri 



n=l 



n=l 



Consequently, we find from (J3.14 ) and ( p. 11 ) that 



K 3 oA 



Sj + iCj 



l-u-i 2Nsm(7rj/N)' 



where 



N 



Sj = J2 B n sin (2n-l)j 



n=l 



7[ 

N 



N 



Cj = J2 B n cos[(2n-l)j 



n=l 



7[ 

N 



ing (p. 16) we obtain 






N C j 
Aj = — arctan — f-, 

2ir Sj 


Kj~- 


_ (s] + qy/ 2 

2Nsm(irj/N) 



(3.12) 

(3.13) 
(3.14) 



(3.15) 



(3.16) 



(3.17) 



(3.18) 



To have Aj = 0, we must choose 6 and such that Cj = 0. From (|3.17 ) and 
( |3~T3|) , we find if 6 + <p = -n/N, then B n = -B N _ n+1 and Cj, Aj = 0. 



3.3. Alternative Approach 

Even though the results in the previous subsection suffice, for future reference 
it may be useful to provide more explicit formulae using a normalization first 
used by Baxter.^ 
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If the average energy for a bond is to be zero, the natural normalization 



is 



N-l 



Yl W(n) = 1. 



(3.19) 



n=0 



We may then rewrite the W(n) in (|2.13| ) in terms of the variables 6 and 
given in (|3.3|) explicitly as 



N-l 



W(n) = J] f m a — , 



where 



m=0 



, _ S^ + g ) _ , R 



sin ffl + 



A/ 



«m,n = | sign(n + | - m) 



N-l 

logW(n) = J2 a »v™ log/„ 

m=0 



(n + \ — m) 
N 

g(n) 



also using the identity 

sin(iVT) = 2 iV - 1 JQ sin(# + 



7V-1 



n=0 



727T 



jv-i 

n=0 



np 



CJ 



?np 



1 -up 



, if P 7^0, 



and its inverse is 



— l uj- np u mp 



E 



{ N 1-ojp 



a r 



We can then rewrite 



N-l 



N-l 



£{n) = -k B T \ogW(n) = --^ a m , n log/ m = JI ^w 

^ m=0 



In 



1=1 



(3.20) 

(3.21) 

(3.22) 
(3.23) 

(3.24) 



It is easy to show, replacing u by x and then taking the limit x — > u>, that 

0, if p = 0, 



(3.25) 



(3.26) 



(3.27) 
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with 

F„ = 

N 



i N-l 

Ei = -J2 £ ( n )^~ ln > E o = ^ (3.28) 



n=0 

see also ( |2.2| ), together with the inverse results 



1 N ~ 1 uj~ ml mi 21 ' 1 *)/ 4 N ~! 

^-p„?,r^ bE/ --^i s,r"' iosU ' (3 - 29) 

-Sjv— z — ^ ) Em = E = 0. (3.30) 



These results also agree with ( |3.16|) of the previous subsection. 



3.4. Superwetting in the Ground State 

Let us say that the interface is superwet at T = 0, if we can find a suitable 
relabeling of the state differences n such that 



£{n) - 8(0) = n(£(l) - £(0)J > 0, n = 0, . . . ,N- 1. (3.31) 

This is equivalent to the second-order difference equation 

£(m+l)-2£(m)+£{m-l) = 0, m = 1, . . ., N- 2. (3.32) 

Imposing the condition, see also (|3.19|) , 



N-l 

£ £(n) = 0, (3.33) 

n=0 

we can solve these conditions, in terms of one constant C, by 

£{n) = -2NC a , n = -C(N - 2n - 1), C > 0, (3.34) 

together with its Fourier coefficients 

U {2l-N)/A 

E '=- C iWS)= £ » l = i,...,w-i, (3.35) 

where we have used (|3.22|) and (|3.25|) . We note that the results ( |3.34| ) 
and (|3.35|) correspond also to the zero-temperature superintegrable© chiral 



Potts model and its corresponding quantum chain. & 
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4. INTEGRABLE THREE-STATE CHIRAL CLOCK MODEL 

In this section, we take a closer look at the 3-state integrable chiral clock 
models. For N = 3, we use the chiral clock representation fl2.12|) to write the 
Boltzmann weights as 



W(n) = exp 
W{n) = exp 



, 27T 

2i<Ticos(— (n + Ai) 

2n 

2i^icos ( —(n + Ax] 



Letting 



(4.1) 
(4.2) 



Wn = W(-n)/W(0), w n = W(-n)/W(0). 
and writing K = 2K\ and K = 2Ki, we find from (|4.1| ), 

w 2 /wi = exp[-v / 3^sin(27rA/3)], w 2 w x = exp[-3ATcos(27rA/3)], (4.3) 

where we have dropped the subscript 1 in A. Consequently, we find 

1 



A' = 2A'j - -^\og\w 1 w 2 ) + 3\og\w 2 /w 1 ) 



and 



A = Ai 

From Q2.21|), we find 



A arctan ( ^QS(w 2 /w 1 ) \ 

cLICLcLIl I " - - I . 

2tt \ log(wiW 2 ) J 



(4.4) 



(4.5) 



Wi = Wi[ 



w 2 = w 2 ( 



sin 2 (0) sin(# + W) sin( 



§*) 



sin 2 (#) sin(0 + ~7r) sin(0 + |7r) 
sin(0) sin 2 (6 l + \n) sin(0 + hi) 
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1/3 



1/3 



(4.6) 



sin((9) sin 2 (0 + \ii) sin(0 + § 
Similar relations hold for K, A, W\, and W2 in terms of 9 and 0. 

4.1. Ferromagnetic regions 

We want to determine for which values of 9 and (f> the (relative) Boltzmann 
weights are positive real and in the ferromagnetic regime. This means that 
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we have to look for a suitable fundamental domain as the relative Boltzmann 
weights have a few simple symmetries: They are periodic modulo tc in 9 and 



Wi(9 + 7T, 0) = W\{Q, + n) = Wi(0, 0), 

w 2 {9 + 7r, 0) = w 2 (9, 4> + ir) = w 2 {8, 0), 
they invert under the interchange of 9 and 0, 

u>x(0, 9) Wl (9, 0) = 1, w 2 ((f>, 9) w 2 (9, 0) = 1, 
and there is a transformation interchanging W\ and W2, 



(4.7) 



(4.8) 



Wi(§7r-&§7r-0)=<u;2(M), ^(f^r - </», |tt - 0) = wi(0,0). (4.9) 

After some work, one finds that modulo it the ferromagnetic regimes are 
given by 



< w 2 < W\ < 1 



for 



§ < W\ < w 2 < I 



for 



-\^<9<(j)< -\<k-8, 
or 7T - 9 < <p < 9 < |tt, 

-ivr-0< 9 < <f)<0, 



or < <p < 9 < ivr 



(4.10) 



(4.11) 



Here, we can ignore the latter two choices ( |4.11| ) as they reduce to the first 
two choices ( 4.1(J| ) under the interchange transformation (O). 

For the integrable 3-state chiral Potts model, we also need to require the 
integrability condition, see (|2.22| ), 



+ = 6» + 6» + ivr, (modulo ix). 



(4.12) 



Combining this with ( |4.10| ) above, we arrive at a unique fundamental domain 
-\ti<9<(J)<-\ti-9, -|tt < < < -|tt - 0, (4.13) 

where 
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< w 2 = w 2 ( 
< w 2 = w 2 ( 



< wi = wi(9, <p) < 1, 

<wx = wi(e^) < i. 



(4.14) 



Note that the integrability condition (|4.12| ) is fully compatible with the 
transformation (|4.9| ), so that other orderings of W\ and w 2 , and of w\ and 
w 2 , correspond to equivalent domains. 



4.2. Critical region 

We now consider the case with 

K = K, A = pA. 



(4.15) 



Substituting ( |4.6| ) into (|4.4| ) and ([4.5[), and using similar equations for the 
barred variables, one can rewrite (|4.15|) to give two conditions between the 9, 
4>, 6, and 0. As we also have ( |4.12| ), we find there is only one free parameter 
left. 

At the Potts critical point, we find 



K c = K c = \ log(l + V3) = 1.49245929 ■ • • 



A = A = 0. 



(4.16) 



From ( |4.5|) and (4^), we see that for A = pA = to hold, we need to have 



W\ = w 2 , wi = w 2 , or 9 
Now we use ([D| and ( |4.17|) to obtain 



rvr, 9 + 



tTT. 



— logu>i 



logwi 



7TVT - 3 



Kr. 



2 "c or 



-12 71 "' 



9 = 9 



:7T. 



(4.17) 



(4.18) 



Near the critical point, we can make the changes of variables 



^7T - 5(f), 



-V2 71 



9=-±TT 

4 

9 = -]tt 
4 



SO = 5<j) + 5(f) - 59. 



(4.19) 



Substituting these expansions into (|4.6|) , and then using (|4.4| ), ( |4.5|) , and 



Q4TT5D, we find 
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59 



50+(p 2 -l)r c 50 2 + (p 2 -l) 2 r- 



2 ~ 2 <50 3 



+ jk(p 2 - 1) \(l5p A - p 2 (U2 + 64 V3) + 15 



12 



+ (p 2 (250 + 122 V3) - 174 - 78 V$) r 2 
+ (p 2 (110 + 78 VS) + 470 + 222 Vs) r c 
- {p 2 + 1)(261 + 171 V3)] 5^ + 0(5<f) 



(4.20) 



and 



p50 - \{p + 1)(<J0 - 50) = 59 - p6<f> - \{p - 1){60 - 50) 
= ~\p (P 2 ~ 1) ((14 + 8 V3) r\ - (53 + 25 V3) r c + 45 + 18 V3 



x 5<P 3 (l + |(p 2 -l)r c 50)+O(50 5 ), 
where we have defined 

rc =l(l + V3)(l + 1 M^) =1.535044409-., 



(4.21) 



(4.22) 



in order to simplify the first three orders in ( |4.20|) . One can also easily verify 
from ( J4.2U| ) and ( [4.21| ) that 59 expressed as a function of 5<p is of the form 
( |4.20| ) with p replaced by 1/p. 

Substituting ([O]) and ( ggg ) back into flOj) and using 



K -K r 



K c t 
l + t 



-K c t, where t 



T-T r 



we find from ( O ) and ( 4.5|) that 

K c t= -\{p 2 + 1)(3 + V3) r c <50 2 (l + (p 2 - 1) r c 50) + O(<50 4 
and 

TT/ 3- v^ 



A = A/p 



TlK c 



(l + I(p 2 -l)r c ^) + O(,50 3 



(4.23) 



(4.24) 



(4.25) 
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Eliminating 5<p from the above equations we find, near the Potts critical point 
on the integrable curve, where K = 2K\ fa K c (1 — t) = | log(l + \/3)(l — t), 
the result 

t = -C t (l + p 2 ) A 2 + 0(£ 2 ,£A 2 , A 4 ), 

C t = i 7T 2 (2 + (7^3 + 12)K C ) = 9.959616982 • • • , (4.26) 



Since the modulus k is given by (|2.25|) , we can use ( |4.19| ) to ( |4.25| ) to express 
it in terms t = T/T c — 1 as 

k 2 = -C k t + 0(t 3 ), C k = 54 ( 2 + v3)^c = 4.98105Q242 ■■■ , (4.27) 
V ; 2 + {7V3 + 12)K c V ; 

which is independent of p at this order. 

The above expansions for the critical region include the symmetric case 
for p = 1 and the Ostlund-Huse model for p = 0. We do not have such 
universal formulae in the low-temperature regime, which we shall discuss 
next. 

4.3. Low- Temperature Limit 

Considering the three-state chiral clock model with K = K and A = pA 
in the low-temperature limit, we find that the limits T — > and p —>■ 
(or p — ► 1) do not commute. The limit p — > does not reproduce the 
Ostlund-Huse model (p = 0) results, nor does the limit p — > 1 reproduce the 
symmetric model (p = 1) results. We shall have to consider several regimes 
separately. 

First, for p^O and p ^ 1, we let 

^ = -^ + (50, = --7r + <J0, = -#, = --7TH-<$, (4.28) 
and we find the leading low-temperature expansion results 

59fa54>5^ = 54>(^5(t)\ , (4.29) 

5^^5(j>, 59 fx 5<p (^ 5<p\ , (4.30) 
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with 

t = ^(i + rtx) V¥T¥TT= ws (431) 

sm(|(l-p)7r) sm(g(l-p)7r) 

Then we see from (|4.4j) and (f4.6|) that 

K = K^ = log <5?/> -> oo, as 50 -> 0, (4.32) 

3v3 



whereas from fl4.5| ) and ( [4.6[ ) we find that A = pA and 



Further terms to the above asymptotic expansions follow by expanding 59, 
5(j), and 59 in powers of 50, 5ip ~ 50^ fe +fe+1 5 50 fe_1 , and 1/ log £■?/>. This can 
be done for 1 < b < oo, or < p < 1. For 6 — > 1 or p — ► 0, 50 6_1 —>■ 1 and is 
no longer small, so that different asymptotic expansions arise. For b — > oo or 
p — > 1, we still have — > ^7r and — > 0, so that there is still another regime 
at p ~ 1 between the above regime and the symmetric case p = 1. 
For the three-state Ostlund-Huse case we have p = A = 0, so that 



sin(3fl) 
sin(30) 



? + ?=-!*, SS = 1. (4.34) 



Using also (|4.13| ), we can then solve 6 1 and in terms of 9 and 0, namely 



9 = -fK-\{9-4>), <j> = \{9-<f>). (4.35) 

In the low-temperature limit, with and 9 near — ^7r, 

0= -±77 + 50, 9 = -\u + 59, (4.36) 

we can find in a systematic way, from K = K, the asymptotic expansion 
for 59 in terms of powers of 50, 50 , and 1/ log 50, in agreement with what 
one would expect from ( |4.30|) and (|4.31| ) with 6 — 1 = 0. The asymptotic 
expansion so obtained can be equivalently (and more economically) expressed 
as 
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log 59 = \og{5(j)5i)) + ]-5(p- y/3 6xjj 



+ 



'4 + v/3 



Zi L», 



1 + ^*0+^V 



+ 77,+ 



2 T r 2 

'5 + 7^ 2 + v 7 ^ v 7 ^ 



18 ' 4L, 2 



6W3 a , a7 o 27-29^, ■■:• 



2 lip 



21 



8cj) 2 8ip 



rp 



where 



lj/j = log 5i/j = V3 log (5(j)/\/3), 



5^j = (5(f)/V3 



v/3 



Therefore the leading order is 

59 « 50 <tyi 



a+%/3 



(4.37) 



(4.38) 



(4.39) 



which is quoted below Q2.66| ) in subsection [2"1)| . Now ( |4.4j ) to (|4.6| ) can be 
used to give 



A' 



|V3log# + ---, A 



1 
4 



Ayr 



e -3X/2 + 



(4.40) 



which was plotted as line 3 in Fig. 7b. 

For the symmetric case with p = 1, the low-temperature limit is given by 



■¥ 



9 = 9 
with 59 ^ 0. We find from 



6" + Itt = -Jtt 

6 6 



to 



A = -§V3log<W + --- 



<J0, (4.41) 

, the asymptotic expansion 



which was plotted as line 1 in Fig. 7b. Thus A decreases from ~ linearly as 
T increases from zero, which is a behavior very different from that shown in 
( |4.42 ) for the Ostlund-Huse case, or from ( 4.33 ). 

The symmetric case can be extended by replacing the products in ( 2.54J) , 
in leading order, by their limiting values and assuming that the front factors, 
which are powers of 1/A and 1/A, are small. The weights then take on the 
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form ( |2.59| ), which was also considered by Baxter. We shall call this case the 



"quasi-symmetric" case. More precisely, with 



9 = -±ir + 89, 6 = -±tc + 56, <f>- 


> 0o ¥" or 


-gff, 


(4.43) 


we find from (|4.13|) that 








= -1^-0 + 59 + 59- 


► -iff -00 




(4.44) 


Now, if we let 








sin 2 ( iff - ( 

n — " 


») 




(4.45) 


sin(g7r)sin(g7T + O ) 


sin(-0 o )' 


and 








sin(-0 o ) 

s = — log -, , s = — 

sin(g7r - O ) 


sin(^ff 
log— -f— 


+ 0o) 

-00)' 


(4.46) 



then we find from ( [4.6|) for (5# and <5# small, the simple formulae 

w n « e" s (# ^) n/3 , W n « e~ 3 (# <0)™ /3 , for n = 1 or 2. (4.47) 

Here we have ignored terms that are one order higher in 56 or 56. If we 
substitute ( fi~4TD in (fill) and ((475]), we find 



K«i Mu + 2s) 2 + \u 2 



3V l ; 3 



3 tan(^ff) m 



Ak — arctan 



6 



K~\\l{u + 2nY + \u\ 








(4.48) 


As 


3 tan(^ff)ii 

i — arctan , 

2tt u + 2s ' 


(4.49) 



2tt u + 2s 

where 

M = -log(^^), u = -log{gS9). (4.50) 

Formulae ([4.48|) and (|4.49|) are correct up to exponentially small corrections 



in the temperature T, as we have ignored terms of order 56 and 56. Therefore, 
u = | V27K 2 - 3s 2 - |s + exponentially small, (4.51) 

and a similar equation with w, K, and s replaced by u, K, and s. Expanding 
these to a few orders in 1/K and 1/if, we find 

- = I^K-ls-^ 2 V3^ + 0^y (4.52) 
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and a similar equation with bars for u. Substituting ( |4.52| ) into ( |4.49| ) , we 
find 

1 s s 3 ^ f s 5 \ 

A = -- — — , + 



4 2nK lOSirK 3 {k 5 )' 

1 5 S 3 ( S 5 ' 

4 2^ \08nK 3 + YE' 



A^-^-^-I + Oi (4-53) 



These again exhibit linear behavior in T, as in the symmetric case. In fact, 



for (j) = (j) = -^tc and 8 = 8, we have s = s = log 2, so ( |4.53| ) reduces smoothly 



to flP3) . 



On the other hand, when — > —-^n, <p — > 0, we have s — > and s — > oo. 

o 

Then A tends to j extremely fast, whereas A tends to j very slowly. In the 
limit, there is a crossover from the linear behavior of ( |4.53| ) to the exponential 
behavior of the first low-temperature regime described by ( f4.28| ) to ( |4.33| ). 

Finally, solving <fio from s = 7 s and 8 from u = 7 u, we find from ( |4.48| ) 
and ( |4.49| ) that K = 7 K and A = A to all algebraic orders in T. So 
the quasi-symmetric case is closely related to integrable 3-state chiral clock 
models with symmetric chiral field, A = A, and we can use this as a 
definition, suitable also outside the low-temperature regime. 

5. SYMMETRIC iV-STATE CHIRAL POTTS MODEL 

In this section we study the symmetric case for general N and give graphs 
for the chiral clock model parameters. 

5.1. Symmetric Lattice 

For the symmetric case with W = W, the integrability condition ( [2.181 ) or 
( |2.22|) becomes 

<P = 8 + —. (5.1) 

V 2N K J 

Substituting this equation into Q3.13J ), we find from (|3.17|) and (|3.18|) that Kj 
and Aj depend on only one parameter 8. We only need to study the regime 
(fundamental domain) with 

K „ 37T . 

< 8 < . 5.2 

N ~ ~ AN K ' 
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First, we can see from ( |3.13| ) that if we let 9 — > 9 + tc/N, then B n — > B n+ \. 
Consequently, combining the first part of ( |3.16| ) and the second equation of 
( |3.15p , we find that KjU Aj — >■ uj^KjUJ K This means that by shifting the 
domain to < 9 < hn/N, the amplitudes Kj remain unchanged but Aj — > 
Aj + j. Secondly, we note the reflection symmetry 9^—9 — \9/N, changing 
B n — > —B_ n . Thus we can shift the domain to —jir/N > 9 > —^tc/N, with 
the amplitudes Kj remaining unchanged but now Aj — > — Aj. Finally, for 
— ^n/N < 9 < (modulo n/N) the Boltzmann weights are not all positive, 
as can be seen easily from (|3.11|) or (|3.20|) , with condition (|5.1| ). 



For 9 I —n/N, we can see from ( |3 . 1 3[ ) that B\ — > — oo, hence Kj —> oo 



and temperature T — > 0. From (|3.15| ) and (|3.16|) we find 



Since the ratios are finite, we obtain the zero-temperature results 

9A i i / n Ka sin(7r/A0 

- 2A ' = -- J . A, = I(JV-2A ^ = -^j- (5-4) 

Comparing with (|3.35|) we see that we are at a superwetting point. 
For the other end of the fundamental domain, 

3-7T 

9 = 0* = -^, (5.5) 



we find using (|3.13| ) that B n = —B N _ n+1 . Hence, we conclude from ( |3.17| ) 
that Cj = 0, or Aj = 0. In fact, by putting 9 = 6*fz — A, it is easy to show 



that B n (X) — >■ — i?7v-n+i(— A) for A — > —A. Consequently, from (|3.17|) we find 



that the Sj are unchanged but Cj — > —Cj, as A — > —A. This means that the 
Kj are unchanged, but the Aj flip signs as A changes its sign. 

To summarize, we find, as we increase 9 from —n/N, that the integrable 
line plotted as a function of 1/Kj versus Aj starts from the zero-temperature 
values given in ( |5.4|) and the Aj decrease and the 1/Kj increase and for 
9 = 9-pz, we have Aj = and the 1/Kj reach their maximum values, 

1 [1 * N] f ?tt\ sin((n-i)7r/iV) 



Nsm(7TJ/N) ^ V v y iVy""° s in((n-f)7r 



i' 
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This corresponds to the symmetric case of the Fateev-Zamolodchikov self- 
dual solution, and is therefore critical. Needless to say that (|5.6Q gives the 
critical temperature for the symmetric solvable chiral Potts model. As 9 
increases further the values of Aj are now negative. The curves are symmetric 
with respect to the vertical axis. 



5.2. Graphs for Integrable Symmetric Case with N=4, . . . , 7 



VKl,2 




Fig. 8a. The 4-state chiral Potts model: \jK\ and I/K2 versus Ai. 



For iV = 4, using (|2.12| ), we have 



^ = 2K X cos Q(jn + A 1 ))+ K 2 (-1)". 



(5.7) 



In Fig. 8a we plot \/K\ and I/K2 versus Ai and in Fig. 8b we plot 1/K 2 
versus l/K\, which is almost a straight line with slope approximately given 
by \f2 ~ 1.41. At temperature T = 0, we have from (|5.4j) that A\ = | and 
K2/K1 = |v2- Therefore, if we let k-oTKi = 1 in the limit T — > by a 
suitable choice of units, then the energies (|5.7| ) for the different states are 



5(0) 

£(2) 



IV2, 



IV2, 



£{l) 

8(3) 



-IV2. 



(5.8) 
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2.5 3 y Kl 



Fig. 8b. N = 4: l/i^ 2 versus l/Ki. 



Hence, 8{—n) — 5(0) = n\/2. This means that at zero temperature the 
integrable line is at a superwetting transition. As T or 1/Kj increases, Ai 
decreases. When Ai = 0, we find from (|5.6|) that 



K lc = 0.3029227993, K 2c = 0.2203433968, K 2c /K 



lc 



0.7273912604. 
(5.9) 
Comparing with the ratio K 2 jK\ = \\pl at zero temperature, we know that 
the curve in Fig. 8b is not a straight line. We may magnify the effect by 
plotting in Fig. 8c the ratio K 2 jK\ versus Ai and indeed we find it not to 
be constant, with k jc <0.02 for Kj defined in ( |2.28| ). 

For different N, we find more or less the same situation as for N = 4. 
For N odd, there are equal numbers of Kj and Aj with 1 < j < |(iV - 1). 
At T = 0, we find then from (WM that Ai 

Aj, not counting Aijy = 0. Now at T 



_(jv_i i - i i s ^ ne smallest of the 
For even N, there are ^A^ different amplitudes Kj and \N — 1 angles 

0, Ai Ar _ 1 = \ is the smallest. 

Then at T = 0, 



Again without loss of generality we can let k^TK\ 
we have 

^ sin(7r/A0 

3=1 



S(n) 



sm(irj/N) 



sin 



(2n 



' N 



(5.10) 
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Fig. 8c. N = 4: K 2 /K 1 versus Ai. 
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Fig. 9a. The 5-state chiral Potts model: \jK\ and I/K2 versus A2. 
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Fig. 10a. The 6-state chiral Potts model: 1/K\, I/K2, and I/K3 versus A2. 
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Fig. 11a. The 7-state chiral Potts model: 1/K\, I/K2, and I/K3 versus A3 
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From this it is straightforward to show that 

£{-n) - 8(0) = 2nsin(vr/iV). 



(5.11) 



It also follows by comparing ( |5.11| ) with ( |3.34| ) and ( |5.1(J| ) with ( |3.35| ). Hence, 



it is at the endpoint of a superwetting line. In Fig. 9a, (respectively 10a, or 
11a), we plot l/Kj with 1 < j < [^N] versus the smallest angle Ari( Ar _ 1 xi for 
N = 5, (6, or 7). In Fig. 9b, (10b, or lib), we plot l/Kj with 2<j< [|iV] 
versus 1/Ki for N = 5, (6, or 7). Again we find that the curves look very 
much like straight lines. 




i«i 



Fig. 9b. N = 5: 1/K 2 versus 1/K X . 

For 9 = 9-pzi we have Aj = 0, and the Kj given by ( |5.6| ) have the numerical 
values 



A^= 7 



K lc = 0.2878960239, K 2c = 0.1845136297, 

K lc = 0.2793167217, K 2c = 0.1675087565, K 3c = 0.1468955978, 
K lc = 0.2738160609, K 2c = 0.1577390412, K 3c = 0.1288523620, 

(5.12) 



which are the critical-temperature values. The ratios are given by 
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Fig. 10b. N = 6: 1/K 2 and 1/K 3 versus 1/i^i. 
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Fig. lib. N = 7: 1/K 2 and 1/K 3 versus 1/i^i. 



55 



N 
N 
N 



K 2c /K lc = 0.6409037097, 

K 2c /K lc = 0.5997090164, K 3c /K lc = 0.5259105038, 

K 2c /K lc = 0.5760766578, K 3c /K lc = 0.4705800001. (5.13) 



Comparing these numbers with the ratios at zero temperature, 



TV 

N 
N 



Kz/Ki = 0.6180339887, 

K 2 /K x = 0.5773502692, K 3 /K 1 = 0.5, 

K 2 jK x = 0.5549581321, K z jK x = 0.4450418679, (5.14) 



we find that the curves are not really straight lines. 




Fig. 9c. N — 5: A 2 versus Ai. 

In Fig. 9c, (10c, or lie), the A.,- for N = 5, (6, or 7) with 2 < j < [±N] 
are plotted versus Ai. These curves also look like straight lines. In Fig. 9d 
we have plotted the ratio K 2 /Ki versus A 2 , while in Fig. 9e, we have plotted 
the ratio A1/A2 versus A 2 , both for iV = 5. Again we find the ratios not to 
be exactly constant, with Kj, Sj <0.02. 
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Fig. 10c. N = 6: A 2 versus Ai. 




Fig. lie. N = 7: A 2 and A 3 versus Ai. 
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Fig. 9d. TV = 5: K 2 /K 1 versus A 
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Fig. 9e. N = 5: Ai/A 2 versus A 2 . 



6. LOW T BETHE ANSATZ FOR DIAGONAL INTERFACES 

In this section we shall present our Bethe Ansatz calculations with free 
boundary conditions&Hl for interfacial tensions of diagonal interfaces in the 
low-temperature limit, using an extension of the Bethe AnsatzO similar 
to the ones introduced by Gaudin.®. These calculations also confirm our 
dilogarithmic integral conjecture fl2.76|) for the entropic corrections (linear in 



the temperature T) to the ground-state interfacial tension results. 

Using Baxter's notation,®® we start with the two coupled eigenvalue 
equations 

T g -y = T g x, f g .x = f,y, (6.1) 

where T q and T q are the two different diagonal transfer matrices. We solve 
this pair of equations, by assuming 

3?m,...,n r 9[ n li ■ ■ ■ i n r\{,li ■ • ■ i SrJ) 

1 < n% < ri2 < • • • < n r < L, 

y ni ,...,n r = g{ni + i,...,n r + ||Ci,...,Cr), 

1 < m < n 2 < • • • < n r < L - 1, (6.2) 

with a common function g which is a linear combination of exponentials. Here 
rii, . . . ,n r stand for the positions of the interfaces in a given diagonal. We 
exclude overhangs, as in the Muller-Hartmann-Zittartz approximation,® 
which is asymptotically exact at low temperatures. If a sequence of, say, 
m positions coincide, or n 3 ■, = . . . — ra J+m _i, then the interfaces merge to a 
single m-step interface. The (\,...,(r are complex parameters, also often 
written as Q = exp(ikj). In terms of these, the eigenvalues are 

T q = f q = (w.w^ 2 n (C7 1/2 + Cf). (6.3) 



12 Choosing fixed boundary conditions for the spins in the first and last rows, corresponds 
to free boundary conditions for the domain walls which live on the dual lattice, which 
is equivalent to demanding that a domain wall can not move outside the system as is 



expressed in (|6.8|) and (3.9) below 
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From these, we can calculate the interfacial tension per two bonds via 

e r = -k B T log(T,f,), (6.4) 

in agreement with Baxter's normalization. 

For the function g we choose the Bethe Ansatz form 

g(mx, ...,m r |Ci,---,Cr) 



£j/2 nn mc ■"'"' 

l<j<k<r 

e nn ^c;?,c^) 1/2 ric;r j ', ( 6 - 5 ) 



* e ••• e n(-cr £i/2 nn ^ccr'.c* 

ei=±l e r =±l j=l l<j<k<r 



X 

P6S r l<j<£:<r j=l 

with the usual two-body scattering function 

^- ^TOt| ^ S 2„A S Q-l)-'. (6,, 

This form guarantees that the "two-interface scattering" is consistent, or 
(w2 — wl) g(nj^n,rij + i^n) — WiWig(nj—>n+l,nj + i-^n) 
+ (W 2 - wl) g(rij^n+ l,nj+i—m+ 1) = 0, 

for j = 1, . . . r — 1. (6.7) 

One can next verify that three-interface and higher scattering processes are 
consistent because of the special form ( p.59|) , which relates to weights of 
fusion models generated from the six- vertex model.® This is also true for 
the case with A ^ A, for which ( |6.5|) need be amended with trivial A/ A 
powers which break the left-right symmetry. 

In (6.5), the last line corresponds to the usual Bethe Ansatz,0 whereas 
the summations over the e's correspond to all waves reflected at the two 
boundaries. The coefficients are chosen such that the first "left" boundary 
condition, 

#(mi = i,...,m r |G,...,Cr) =0, (6.8) 

is automatically satisfied. 

In order to solve the other "right" boundary condition 

5f(mi,...,m r . = L + ||Ci,--.,Cr) = 0, (6.9) 
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we have to require the "Bethe Ansatz equations" 

Ck L U{ A «d,OdM<j\Ck)) = l, k = l,...,r. (6.10) 

We need to solve these equations ( p.l(J| ) in the thermodynamic limit L — > oo 
at finite r and substitute the solutions in ( |6.3|) and ( |6.4|) . 
The large- L solution giving the largest eigenvalue is 

cob((|-[^])A) . 

= .»(«-[¥])aV ' '-" ( ' 



with [x] the integer part of x. For even r (|6.11| ) reduces to 



cos Us- j + l)A 

r = 2s, Ciy-i = C« = — 4; rryA j = l,...,s, (6.12) 

cos(Js-j)AJ 

whereas for odd r it reduces to 

cosUs-j + |)A) 

r = 2s + 1, C2J-1 = C 2j = j- : — f , j = 1, . . . , s, 

cos^s-j + 2)AJ 

C2.+1 = 1. (6.13) 

From these results we immediately obtain fl2.75|) with a r given by the last 
member of Q2.76J ). We see that we are precisely at the boundary where the ( 
pairs change from complex conjugate pairs to unequal real pairs, which may 
be related to some pre-wetting phenomena on the integrable line. 

7. NEW MODEL 

There is already a great deal of exact results obtained for the chiral Potts 
models. Yet not much of this has been utilized. One of the reasons is that 
the temperature dependence in the integrable chiral Potts model, as given 
by ( |2.7|) through (2J3) is not very transparent. It is rather vague what the 



critical temperature T c is, thus making it harder to use the exactly known 
exponents and other results for N > A. 
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Moreover, one would like to have a model which includes the chiral Potts 
model as a special case, yet not so general as (|2.2|) which has too many 
variables. Such a model, with N > 3, can be used to describe some 
interesting physical phenomena, such as commensurate-incommensurate 
phase transitions and wetting phenomena. It may not be exactly solvable 
in most regions. Hence, numerical studies would have to be done on it. 
However, the model may be in the same universality class as the integrable 
chiral Potts model, and may be more interesting for N > 4 than the chiral 
clock model which has only a single cosine term in (|2.12 ). 



7.1. New 2-Parameter Chiral Potts Model for General N 

We therefore propose a new model with nearest-neighbor pair interactions 
given by 

This has only two variables E and A. For A = |, it is identical with the 
pair-interaction energy ( |5.10|) of the integrable chiral Potts model at zero 



temperature and is thus on the superwetting line there. For A ~ 0, we find 
£(0) is the ground state. In fact, it is easy to show that for A < A CI , with 

A - 1 ( N + 1 ) (72 ) 

and at low temperatures the system is in the ferromagnetic phase (or W(0) 
is the maximum). Within this regime, for A < |, the system is unwetted, 
namely e p + e r _ p > e r ; at A = |, it goes through a superwetting transition, 
with e p + e r _ p = e r ; while for A > |, the system is wet, with e p + e r _ p < e-. 
At A = A CI , we find that £(0) = £(-1) as in the chiral clock model.©© 
Therefore, the ground state is highly degenerate and we believe that, at 
this point, the system goes through a commensurate-incommensurate phase 
transition. 

7.2. Superwetting and Ground-State Degeneracy for N—5 

To make all these statements clearer, let us take iV = 5 as an illustration. 



We plot, in Fig. 12, the energy £{n) given in (|7.1| ) as a function of x = n/N, 
for different values of A. 
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-£(«) 4 




x=n/5 



Fig. 12a. Energy —£(ri) (or interfacial tension at zero temperature) versus x = n/N = n/5 
for the new model with N — 5 and A = i : ground state in the "dry phase" as no wetting 
can take place. 

In Fig. 12a, we show the situation for A = | < |. We find that the 
increments in energy - - or the interfacial tension at zero temperature - 
e r = £{—r) — 8(0) are given by 



e 4 
(2 



4.574329190, 
3.315920618, 



£3 



4.209056927, 
2.036147842 



(7.3) 



and are shown in Fig. 12a as the lengths of vertical straight lines indicating 
the energy differences between the 8(—r) and the ground state energy 8(0). 
By examining the lengths of these straight lines, we find that e p + e r _ p > e r . 
Therefore it is energetically unfavorable to have two interfaces with interfacial 
tensions e p and e r _ p instead of one with e r . Hence, the system is not wetted. 
In Fig. 12b, we plot the energy 8 for A = |, and we find e r = rt\ = 
2r sin n/5. This means that it is energetically neutral to have one interface 
with e r or more than one, even r interfaces with e\. This is what we have 
called superwetting. 



In Fig. 12c, we plot 8 for A = ^ > \. We now find 



(2 



10 

4.655997407, 
1.703676245, 



e 3 = 2.965539502, 
ei = 0.6789144637 



(7.4) 
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x=nl5 



Fig. 12b. Energy — £(ri) for TV = 5 and A = j: the ground-state is at a (super)wctting 
transition. 



-£(«) 




Fig. 12c. Energy — £ (n) for N — 5 and A = y^: ground state in the wet phase. 



64 



and also that re\ < e p + e r _ p < e r . Consequently, it is energetically more 
favorable to have r interfaces with e\ than just one with e r or another number 
less than r. This means the system is wetted, allowing only one kind of 
interface with interfacial tension ei. 



-E{n) L\ 




x-n/5 



Fig. 13. Energy — £ (n) for N 
W(0) = W(N- 1) here. 



5 and A 



the ground state is highly degenerate as 



We plot £ for A = § in Fig. 13, and we find S(N - 1) = S(0). This 
is the same situation as in the chiral-clock model, namely that the ground 
state is highly degenerate for certain values of A. Hence, near A = | and at 
low temperature, there exist floating incommensurate phases. The model 
( |7. 1|) can then be used to describe commensurate-incommensurate phase 
transitions. However, in these regimes, exact solutions do not exist and 
numerical methods have to be be employed. 



7.3. Comparison with Integrable Model for N=5 

Comparing with the values given by ( |5.4j ) for the integrable model, we 
find that the new model given by ( [7. 1|) is integrable at zero temperature at 
the special value A = \. On the other hand, unlike for the new model, 
the ratios KjjK\ and Aj/Ai for the integrable cases are not constants, as 
shown in Fig. 8 to Fig. 11. We believe that the new model may very well not 
be integrable at any other point. On the other hand, as the ratios for the 
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x=nl5 



Fig. 14a. Comparison of the energies — £ of the self-dual integrable chiral Potts model 
and the new model with A — for N = 5; the two curves are almost on top of each other. 



-e(n) A 




x=nl5 



Fig. 14b. Comparison of the energies — £ of the integrable chiral Potts and the new model 
at A = .0977875 for N — 5; again the two curves are almost on top of each other. 



66 



integrable cases are almost constants, it is interesting to compare them. In 
Fig. 14a, we plot for N = 5 the nearest-neighbor interaction as a function of 
x = n/N for the integrable case at the self-dual point with 

Ax = A 2 = 0, K x = 1, K 2 /K x = 0.6409037097, (7.5) 

together with that of the new model at A = 0. Also, in Fig. 14b, we plot the 
8 for the integrable model evaluated at 9 = — 97r/50, with 

Ax = 0.2943916857, A 2 = 0.09778750304, 

K x = 1, K 2 /K x = 0.6404446571, (7.6) 



and the £(n) of evaluated at A = A 2 = 0.09778750304. We find that 
the curves are almost on top of each other. 

As it is well known that by changing the strengths of the interactions, 
without changing the symmetry (here N) or spatial dimension, the exponents 
are universal, it is likely that the integrable model and the new model have 
the same critical exponents for the same N. On the other hand, the models 
must be different in some other sense and it would be interesting to find out 
more about them. 

It may also be interesting to study variations of model ( [7.1[ ), which 
agree with the symmetric integrable model at a non-zero temperature and to 
investigate if and how superwetting sets in. In the most general chiral Potts 
model there are enough parameters to allow the existence of solutions of 
the conditions for a superwetting line. It may be of interest to see, by low- 
temperature expansion for example, if there are solutions with real positive 
Boltzmann weights and what their physical implications are. 
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